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Abstract. The interaction between a viscous fluid and an elastic solid is 
modeled by a system of parabolic and hyperbolic equations, coupled to one 
another along the moving material interface through the continuity of the 
velocity and traction vectors. We prove the existence and uniqueness (locally 
in time) of strong solutions in Sobolev spaces for quasilinear elastodynamics 
coupled to the incompressible Navier-Stokes equations. Unlike our approach 
in [3 for the case of linear elastodynamics, we cannot employ a fixed-point 
argument on the nonlinear system itself, and are instead forced to regularize 
it by a particular parabolic artificial viscosity term. We proceed to show that 
with this specific regularization, we obtain a time interval of existence which 
is independent of the artificial viscosity; together with a priori estimates, we 
identify the global solution (in both phases), as well as the interface motion, 
as a weak limit in srong norms of our sequence of regularized problems. 


1. Introduction 

We establish the existence and uniqueness in Sobolev spaces of strong solutions 
to the unsteady fluid-structure interaction problem consisting of a nonlinear large- 
displacement elastic solid coupled to a viscous incompressible Newtonian fluid. The 
fluid motion is governed by the incompressible Navier-Stokes equations, while the 
solid, which can be either compressible or incompressible, is modeled by the cele¬ 
brated St. Venant-Kirchhoff constitutive law (although our method can be applied 
to more general quasilinear hyperelastic models). 

The first fluid-solid interaction problems solved were for the case of a rigid body 
inside of a viscous flow in a bounded domain (see and the case of 

a rigid body inside of a viscous flow in an infinite domain am m, uni). Later, 
the elastic body was modeled with the restriction of either a finite number of modes 
( 0 ) or with a hyperviscous type law for the solid ( 0 , ( 101 ), essentially by the same 
type of Eulerian global variational methods developed in 0. For the steady-state 
problem, which is elliptic in both phases, m solved the case of solid modeled as 
a St. Venant-Kirchhoff material. In m , an Eulerian approach was used for the 
case in which the solid is a visco-hyperelastic material, which is a regularization of 
a hyperbolic model of solid deformation. 

With the exception of our recent well-posedness result for the case of a linear 
elastic solid in [5], there are no known existence results for fluid-structure inter¬ 
action when the solid is modeled by a standard second-order hyperbolic equation. 
This may be attributed to the difficulties associated with coupling a parabolic PDE 
for the fluid with a hyperbolic PDE for the solid through continuity of the veloc¬ 
ity and traction vectors across the moving material interface. As we explained in 
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an iteration scheme between fluid and solid phases fails to converge due to a 
regularity loss induced by the hyperbolic phase (this divergent behavior has been 
computationally noted as well in El), and so we developed a method comprised 
of the following new ideas: first, a functional framework which scales in a hyper¬ 
bolic fashion for both the fluid and solid phases. This scaling leads to additional 
compatibility conditions in the fluid phase (when compared to the use of the classi¬ 
cal parabolic framework), and is absolutely crucial for obtaining consistent energy 
estimates. Second, we developed a regularity theory founded upon central trace 
estimates for the velocity vector restricted to the interface, rather than traditional 
interior regularity arguments which do not work for our problem. Third, we were 
forced to bypass the use of the frozen (or constant) coefficient basic linear problem, 
which requires estimates on one more time derivative of the pressure function than 
the initial data allows, and created a new method wherein the solution was found 
as a limit of a sequence of penalized problems set in the Lagrangian framework. 
The penalization scheme approximates the divergence-free constraint, whereas the 
Lagrangian framework alleviates the difficulties associated with the lack of a priori 
estimates in the solid phase for the frozen coefficient problem; this method indeed 
differs significantly from the classical methods used in fluid-fluid interface problems 
(see for instance E3, HI). 

The fundamental difficulty in extending our result to the case of nonlinear elas¬ 
ticity is the absence of any method of analysis for quasilinear elastodynamics which 
is compatible with the general scheme of |H|, involving a global Lagrangian vari¬ 
ational formulation and the use of difference quotients to track the regularity of 
interface data. We remind the reader that unlike the analysis of elastostatic mo¬ 
tion, direct inverse function theorem arguments cannot be applied directly to the 
case of quasilinear elastodynamics due to the fact that the perturbation term aris¬ 
ing from the nonlinear operator is not an element of the appropriate function space 
for optimal regularity. Alternatively, one might attempt a fixed-point approach, 
wherein a portion of the nonlinear elasticity operator is viewed as a forcing func¬ 
tion coming from a given velocity v, and then try to solve a linear problem for an 
unknown w. The difficulty in this approach stems from the fact that one has to find 
exact time derivatives of elastic energies for the forcing term associated to the elas¬ 
ticity operator, which is complicated by the inner-product of a term involving J'* v 
and a term involving w. This difficulty is overcome in [5|, by a clever and essential 
use of the Dirichlet boundary condition in order to reformulate the problem in a 
non-standard way. As it turns out, the variety of known methods that have been 
used in the well-studied area of quasilinear elasticity, such as those in and m 
for the Dirichlet boundary condition, or in m and i for the Neumann bound¬ 
ary condition, require a priori knowledge of the boundary data regularity, and are 
hence instrinsically incompatible with fluid-structure interaction analysis (in fact, 
the methods devised for Dirichlet conditions do not work for Neumann conditions 
and vice versa). Indeed, of these various methodologies, only Jj| and m use a 
variational approach, the others employing either semi-group techniques as in the 
early work of m in the full space, or the technical paradifferential calculus in m 
for the two dimensional Neumann case. 

In this paper, we develop a new method for quasilinear elastodynamics, varia¬ 
tional in nature, that is compatible with our method in jjj|; we proceed in two steps. 
First, we add a specific artificial viscosity to the solid phase which regularizes the 
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system, thus converting our hyperbolic PDE into a parabolic one, transforming 
the fluid-structure interaction into a fluid-fluid interface-type problem for which 
existence and uniqueness of solutions is already known on a time interval that a 
priori shrinks to zero as the artificial viscosity n tends to zero. Second, and this 
is where the primary difficulty rests, we prove that our specific choice of parabolic 
smoothing renders the time interval (on which a unique solution exists) indepen¬ 
dent of k; furthermore, our a priori estimates allow us to construct a solution by 
weak convergence in strong norms. We note that the use of higher-order operators 
in the artificial viscosity term, while providing the necessary a priori control of the 
regularity of the moving interface, would not yield ^-independent estimates which 
are essential here. Also, as our parabolic reguralization method is not specialized 
to any particular boundary condition, it thus provides a unified approach to the 
classical problem of quasilinear elastodynamics when the solid is not coupled to a 
fluid. 

We now proceed to the formulation of our problem. Let fl C K 3 denote an 
open, bounded, connected and smooth domain with smooth boundary dfl which 
represents the fluid container in which both the solid and fluid move. Let C fl 

denote the closure of an open and bounded subset representing the solid body at 
each instant of time t £ [0, T] with fi-f(t) := fi/fl s (t) denoting the fluid domain at 
each t £ [0,T], Note that in our analysis Q s (t) is not necessarily connected, which 
allows us to handle the case of several elastic bodies moving in the fluid. 

Remark 1. If a function u is defined on all of SI, we will denote uf = u 1^/ and 

“0 

u s = u 1q« . This allows us to indicate from which phase the traces on 

r(0) :=fi/(0)nfR(0) 

of various discontinuous terms arise, and also to specify functions that are associ¬ 
ated with the fluid and solid phases. 

For each t £ (0, T], we wish to find the location of these domains inside O, the 
divergence-free velocity field u? (f, •) of the fluid, the fluid pressure function p(t, •) 
on lV(f), the fluid Lagrangian volume-preserving configuration ?/(t, •) : = 

17q —> and the elastic Lagrangian configuration field if(t,-) : f2 s (0) = fig —> 

fl s (f) such that f l = rj s (t, fig) where rj{(t,x) = u^(t,r]^(t,x)) and uf 

solves the Navier-Stokes equations in Qf(t): 

u{ + (t/ • V)w/ = divT^ + ff , 
div u-f = 0 , 

with 

= v Def — pi, (1.1) 

and rj s solves the elasticity equations on fl s (0) 

ff = di vT s + f s , 

with T s = ■ | Tr((V? 7 s ) T V? 7 S — I) I +p (( \/r] s ) T \/r] s — I), and where the equations 
are coupled together by the continuity of the normal component of stress along the 
material interface T(t) := fi s (f) flO-f ( t ) expressed in the Lagrangian representation 
on To := T(0) as 

T s N = [T f or/] [(V?/)" 1 N ], 



4 


D. COUTAND AND S. SHKOLLER 


and the continuity of particle displacement fields along To 

f S 

v = n , 

together with the initial conditions u( 0, x ) = Uo(x), and 77 ( 0 , x) = x and the Dirich- 
let (no-slip) condition on the boundary dtt of the container uf = 0, where v > 0 is 
the kinematic viscosity of the fluid, A > 0 and p > 0 denote the Lame constants of 
the elastic material, N is the outward unit normal to Fq and Def u is twice the rate 
of deformation tensor of u, given in coordinates by u l ,j +u J All Latin indices run 
through 1,2,3, the Einstein summation convention is employed, and indices after 
commas denote partial derivatives. 

We now briefly outline the proof. As the solid and fluid phases are naturally 
expressed in the Lagrangian and Eulerian framework, respectively, we begin by 
transforming the fluid phase into Lagrangian coordinates, leading us to the system 
of equations ed and, as in (3J, we work in an hyperbolic framework in order to 
accomodate the dual nature of the problem (parabolic in the fluid and hyperbolic 
in the solid). 

In order to solve iTHl) . in Section Q we first add a particular form of artificial 
viscosity to the quasilinear hyperbolic equation in the solid, transforming the hy¬ 
perbolic phase into a parabolic one; specifically, we add the term where L 

denotes the linearized (about the identity) elasticity operator and ry is the material 
velocity. We hence obtain an interface problem that is parabolic in nature in both 
phases, and can be thought as a fluid-fluid parabolic interface problem for which 
wcll-posedness is classical, save for the fact that both phases are required to scale 
in an hyperbolic fashion. The time interval of existence [0, T K ] for this parabolic 
system a priori shrinks to zero as n —> 0. 

In Section 0 we establish ^-independent estimates on the solutions v K to the 
regularized parabolic problem on the time interval [0, T„] by identifying exact time 
derivatives of elastic energies, and establish regularity of the interface. A direct 
fixed-point approach for (El does not appear to yield these exact time derivatives 
for the elastic energy, whereas the the regularized problem ED does indeed lead 
to them. An essential key for obtaining estimates independent of k inside the solid 
is Lemma ED Whereas the trace estimates could be carried with other choices of 
artificial viscosity, we absolutely need the special choice made in our analysis in 
order to recover the regularity inside the solid independently of k. In particular, a 
different choice of a regularizing operator either of the same order such as —A rj t or 
of higher order such as L 2 rj or L 2 r/t would not provide ^-independent estimates. 

In Section El we then explain how our estimates allow the construction of solu¬ 
tions v K on a time interval independent of k, still with energy estimates independent 
of k. The existence of a solution to ED then follows by weak convergence as n —> 0. 

Uniqueness is established in Section ED in the same functional framework used 
for existence. 

As our method seemingly requires more regularity on the initial data in the solid 
than it should, due to the artificial viscosity in the compatibility conditions, we 
explain in Section IT2l how this extra regularity can be removed, thus leading to the 
result with optimal regularity. 

SectionEHis dedicated to the case where the incompressibility constraint is added 
to the solid. The additional difficulty with respect to the compressible case comes 
from the fact that we control the velocity uniformly in n in function spaces which 
possess less regularity than in the fluid, whereas the pressure is controlled uniformly 
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in the same regularity spaces in both phases. Also, we cannot use Lemma |ti.ll in 
the most optimal form for the regularity of the pressure in the solid phase. 

2. NOTATIONAL SIMPLIFICATION 

Although a fluid with a Neumann (free-slip) boundary condition indeed obeys 
the constitutive law o, we will replace for notational convenience o with 

Tf = v'S/u* — pi; (2-1) 

this amounts to replacing the energy / Def : Defv by / Vw ^ : Vu, which is 

Jn* ' Jn s 0 

not a problem mathematically due to the well-known Korn inequality. Henceforth, 
we shall take O as the fluid constitutive law. 

3. LAGRANGIAN FORMULATION OF THE PROBLEM 

In regards to the forcing functions, we shall use the convention of denoting both 
the fluid forcing // and the solid forcing f s by the same letter /. Since // has to 
be defined in fI (because of the composition with rj ), and f s must be defined in fig, 
we will assume that the forcing / is defined over the entire domain fl. 

Let 

a(x) = [CofV? 7 ^(ai)] T , (3.1) 

where (V?y^( x))j = d^) 1 / dx ^( x) denotes the matrix of partial derivatives of r/A 
Clearly, the matrix a depends on r] and we shall sometimes use the notation a* (ry) 
to denote the formula O- 

Let v = u o rj denote the Lagrangian or material velocity field, q = p o 77 is the 
Lagrangian pressure function (in the fluid), and F = ff o r/f is the fluid forcing 
function in the material frame. Then, as long as no collisions occur between the 
solids (if there are initially more than one) or between a solid and <9f2, the problem 
can be reformulated as 


r) t = v 

in (0, T) x n, 

(3.2a) 

v\ - v{a\a^v\ k ),j +(a k q),k = F l 

in (0,T) x Hq , 

(3.2b) 

a k v\ k = 0 

in (0,T) x Hq , 

(3.2c) 

v t ~ c mjfcz [(?7, m ‘V,j Smj)V,k],l = f 

in (0,T) x fig , 

(3.2d) 

C mjkl (■ r),m -ihj -6mj)v\ k Ni=v v\ k a k a\Nj 

- qa{Nj on (0, T) x T 0 , 

(3.2e) 


a.e. in (0, T) , 

(3.2f) 

O 

$ 

II 

on fl x {t = 0} , 

(3-2g) 

q = Id 

on Q x {t = 0} , 

(3.2h) 


where N denotes the outward-pointing unit normal to To (pointing into the solid 
phase), and 

c ijkl = X6 ij 6 kl + ix(S ik S jl + S il 6 jk ) . 

Throughout the paper, all Greek indices run through 1, 2 and all Latin indices 
run through 1, 2, 3. Note that the continuity of the velocity along the interface is 
satisfied in the sense of traces on To by condition whereas the continuity of 

the normal stress along the interface is represented by E2B. 
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4. Notation and conventions 


We begin by specifying our notation for certain vector and matrix operations. 

We write the Euclidean inner-product between two vectors x and y as x ■ y, 
so that x ■ y = x l y l . 

The transpose of a matrix A will be denoted by A T , i.e., (H T )} = Aj. 

We write the product of a matrix A and a vector b as A b, i.e, (A by = . 

The product of two matrices A and S will be denoted by A-S, i.e., (A-S)j = 

4 s*. 

For T > 0 and k £ N, we set 

Vf(T) = {re £ L 2 (0, T; R 3 )) | d?w £ L 2 (0, T; H k ~ n (Q f 0 - R 3 )), n = 1, 

with V, fe (T) defined with fig replacing 

In order to specify the initial data for the weak formulation, we introduce the 
space 

Ldi V j = {V> £ L 2 (fl; R 3 ) | di yip = 0 in Hq, ijj ■ N = 0 on dtt} , 

which is endowed with the L 2 (fl;M 3 ) scalar product. 

The space of velocities, X T , where the solution to iTOll lives, is defined as the 
following separable Hilbert space: 


X t = {v€ L 2 (0, T; Hq(Q; R 3 )) | (yf, / V s ) £ Vf(T) x V S 4 (T)} , (4.1) 


endowed with its natural Hilbert norm 


v\\~ X - = »v\ 


L 2 (0,T;Hq (Q;R 3 )) 






k 


li 2 (0 ,T;ff 4 -"(Q ( 


|;R3))] • 


n=0 

We also need the space 

Yt = {(v,q) £ X T x L 2 (0,T;ET 3 (^;K))| d?q £ L 2 (0, T; H 3 ~ n (n f 0 -, R))(n = 1,2)}, 
endowed with its natural Hilbert norm 

IK-.9)I| ! v t = IMliv + t, ll3.Nllh(0,r ; „.- (nJi ,), ■ 

n =0 

We shall also need L°°-in-time control of certain norms of the velocity, which 
necessitates the use of the following closed subspace of AY: 

W T ={veX T \ v m £ L°°(0,T;L 2 (H;R 3 )), 




£ L°°(0,T;17 4 -Y^;K 3 ))(n = 0,1, 2,3)} , 


endowed with the following norm 

3 ' r 

IMIwr = IMIx t + IK'«tlll=°(0,T;L 2 (n;R 3 )) + E / ' l ’lli oc> (0,T;ff 4 -"(ng;R 3 )) • 


n —0 
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Finally, we will also make use of the space 

Z T = {(«,g) G W T x L 2 (0,T;iJ 3 (f^;R))| d?q G L 2 (0, T; H 3 ~ n (n f 0 ; R)), (n = 1,2) 

| fc 6L“(0,r;I 2 (ni;R))} , 


endowed with its natural norm 


\\(v,q)\\z T 


2 


INI 


2 


+ EHM 

n—0 


2 

L 2 (0,T;-U 3 - n (fi£;R)) 


+ 11^911 


2 

L°°(0,T;L 2 (Oq ;R)) ' 


Remark 2. Note that our functional framework does not make use of the third time 
derivative of the pressure qttt, even though we do use the third time derivative of 
velocity wttt / this functional framework is necessitated by the fact that the Dirichlet 
boundary condition together with the limited regularity of wttt does not allow us to 
obtain qttt with the appropriate regularity. Note also that we have added the L°°- 
in-time control of q tt in the definition of Zt mostly for a more convenient way to 
prove our theorems, rather than absolute necessity. 

Throughout the paper, we shall use C to denote a generic constant, which may 
possibly depend on the coefficients v, A, p, or on the initial geometry given by 
fl and fig (such as a Sobolev constant or an elliptic constant). For the sake of 
notational convenience, we will also write u(t) for u(t, •). 


5. A FIRST THEOREM 

We now state our first theorem. We impose greater regularity requirements on 
the initial data than is optimal so as to avoid technical difficulties associated with a 
particular type of initial data regularization that would otherwise be necessitated. 
We consider the case of optimal regularity on the initial data in Theorem ll2.ll 

Theorem 5.1. Let fl C R 3 be a bounded domain of class H A , and let fig be an 
open set (with a finite number > 1 of connected components) of class H 4 such that 
fig C fl and such that the distance between two distinct connected components of 
flg (if there are multiple solid components) is greater than zero. Let us denote 
fig = fl D (flg) c . Let v > 0, A > 0, p > 0 be given. Let 

(. fjufttjttt ) G £ 2 (0, T; if 3 (fl; R 3 ) x H 2 ( fl;K 3 ) x H\ f!;R 3 ) x L 2 ( fl;R 3 )), 

(5.1a) 

/( 0) G H\ fl;K 3 ) , f t ( 0) G H\ fl;R 3 ). (5.1b) 


Assume that the initial data satisfies 


u 0 G H\n f 0 ; R 3 ) n H 6 {n s 0 ; R 3 ) n Hq(CI', R 3 ) n L 2 divJ , 
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as well as the compatibility conditions 

[Vitg TV] tan = 0 on To, W\ = 0 = u >2 on 8 f2, vAuq — Vgo = 0 on r 0 , (5.2a) 
[{u[{aUiWXW)Nj)U]t^~ Kao ai t (0)Nj)UWn 

= [c mjkl [(V s ,m -V s ,j -S mj )V S ,k] t (0)m t an 077 T 0 , (5.2b) 

[(v[( a l a l) wf i]tan- [( 2 <?1 tti t (0)iVj+5o a itt(0)-^j')i=i] t a n 

= [c mjfei [(77 s , m -77 s ,i —d m j)v s ,k ] tt (0)TV/] tan 077 Tq, (5.2c) 
vAw{ + v((a : [a't)t(t))u f 0 ,k), j +F t ( 0) - [(K)t(0)g o ),.,- +gi,»]?=i 

= /t(0) + c m -' w [[( 77 s ,m ‘V s ,j —$mj)v s ,k],i ] t (0) 077 T 0 , (5.2d) 

where the time derivatives appearing in these equations and in the following ones 
are computed from any w satisfying w( 0) = uo, 8pw(0) = w n (n = 1,2), and from 
any q satisfying <9"g(0) = q n (n = 0,1, 2), the quantities w n and q n being defined 
as follows. First, go £ H 3 (fl g;R) is defined by 


Ago = div/(0) + (aj) t (0)uQ,j in fl f 0 , (5.3a) 

go = 77 [Vt7q TV] ■ TV on To, (5.3b) 

= /(0) • TV + vAuq ■ N on 8tt, (5.3c) 

and Wi £ 7Jg (Q; R 3 ) n l? 4 (fi§; R 3 ) fl H 4 (Q gj R 3 ) by 

wi = uAuq — Vgo + /(0) in fig (5.4a) 

Wi = /(0) in Qq . (5.4b) 

Note that w\ £ H 4 (flg; R 3 ) since Awi £ H 2 (f2g ; R 3 ) and w\ = 0 on <9f2, w\ = /(0) 
077 To- We also have gi £ H 3 (H gjR) defined by 

Agi = div[z/Awi + F t (0) + K(a* a?M°K,fc )d -(( a i)t(°)9o),j ]?=i] 

+ 2 ( a i)i(°) u, id+( a i)it(°) w od 777 , (5.5a) 

gi = TV • TV + (ai al) t (0)ul\k NjNf\ - g 0 a^(0)TV,TVj 

- c mjkl [( 77 s , m - 77 s ,,- —<5 m j)?7 s ,fe] t (0)Aj • TV 077 r 0) (5.5b) 

^ =-Ft(0) • TV - [(a 4 ) t (0)g o ],j TV £ + t^Awi • TV + z/((afaf) t (0)uo,/ s ),^ TV* 077 911 , 

(5.5c) 


and w 2 € ATq (S7; R 3 ) D H 4 ( fig; R 3 ) D if 2 (f2g; R 3 ) by 

w l 2 = vAw\ + v((aja k ) t (0)uQ,k ),j +Ff( 0) - ((af ) t (0)g 0 ),j -gi,i in , (5.6a) 

W 2 = ft(0) + c mjkl [[(v, m -V,j -$mj)v,k],i] t (0) in£l s 0 , (5.6b) 
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Finally, q 2 G H 1 ^ g;K) is defined by 

Aq 2 = div[(/ o t])t t (0) + v [(af a$w, k ,j ] tt (0) - [((a- )tt{0)q o + 2(aj) t (0)q 1 ), j ]f =1 ] 

+ ^( a i)t(0)w2,j +3(aj)tt(0)w{,j +(aj)ttt(0)uQ,j infl f 0 , (5.7a) 

q 2 = i'[(a k a J l )w fl , k ] tt {0)N j N i - c mjkl [(rf, m -r] s ,j -S mj )r] s lk ] tt (0)Ni ■ N 

- q 0 a i tt (0)NjNi - 2qi a J it (0)NjNi on T 0 , (5.7b) 

=(/ 0 v)tt(0) • N - 2[(af ) t (0)gi],j Ni - [(a^) tt (0)g 0 ],j V + • N 

+ 2v{{a\a k ) t (Q)w\, k ),jN i + v({a\a k )tt(Q)uo, k ),jN i on dfl . (5.7c) 

Then there exists T £ (0 ,T) depending on uq, f, and fig, such that there 
exists a unique solution (v,q) £ Zt of the problem 1.7, S\) . Furthermore, g £ 
C 0 ([0,T];ff 4 (flf;R 3 )nff 4 (£2 s o ;R 3 )nH 1 (fl;R 3 )). 

Remark 3. The remarks appearing in 0 at the end of Section 5 concerning the 
compatibility conditions and forcing functions for the linear elasticity case still hold 
in this setting with the necessary adjustments. In particular, we do not need the 
forcing functions to have the same regularity in both phases. 

6. Preliminary result 
In the remainder of the paper, we set 

L(uY = [c i4kl (u k , l +u l , k )}, j . 

In our limit process as the artificial viscosity tends to zero, we will make use in 
a crucial way of the basic following result: 

Lemma 6.1. Let g £ C°([0, T]; L 2 (fig; R 3 )) and u be such that ut £ L 2 ( 0, T; H 2 (LIq; R 3 )) 
and 

eL(ut) + L{u) = g on [0, T\ x fig. (6.1) 

Then, independently of e > 0, 

IIU M )l|L~(0,T;L 2 (ng;R3)) < ll5l|L“(0,T;L 2 (ng;R3)) + || L(u 0 ) || i 2 (ng;R3) ■ 

Proof. Since L(u) £ C°(0, T; L 2 (Qq; R 3 )), let T' £ [0, T] be such that 
IIU u (D)l|L 2 (ng;R3) = sup ||L(u)|| L 2( n g. R 3). 

[0,T] 

If T' = 0, then the statement of the Lemma is satisfied. Now, let us assume that 
T' £ (0,T]. Let 5 £ (0 ,T') be arbitrary. From 16.1 ll . we infer that 
fT' fT' 

e 2 l|i(wt)llL 2 (ng;R3) + / IIU u )llL 2 (ng;R3) 

JT'—S JT'—S 

t' 

+ 6 [IIU u (T))ll! 2 (f2g;R3) - II L(u{T' - 5))|||2( n «. K 3)] = / ll5lll 2 (ng;R3)- 

J T'—S 

From the definition of T' we then infer that for any <5 £ (0, T'), 

[ II II L 2 (Oq ;R 3 ) — f II5II L 2 (Oq;R 3 ) ’ 

JT'—S JT'-S 
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which after division by S gives at the limit 5 —» 0: 

||I/(-u(T / ))||| j 2( n g ;R 3) < ||sOOIII 2 (q* ; ir 3 )> 

which concludes the proof of the Lemma. □ 

Remark 4. It should be clear that Lemma fO applies to a more general class of 
linear operators than L. 


7. The smoothed problem and its basic linear problem 

As we described in the introduction, we cannot find an appropriate linear prob¬ 
lem whose fixed-point provides a solution to (IQ) . We are thus lead to introduce 
introduce the following (parabolic) regularization of 13.21 1, with the artificial viscos¬ 
ity k > 0 : 

v l t -u(a : l afv l , k ), j +(a^q), k = f 1 on in (0,T)xfl£, 

(7.1a) 

a*V, fe = 0 in (0,T)xOj, 

(7.1b) 


vi-K[c ijkl v k , l ], j +N(n) i = f i + K h i 

in (0 ,T)xflg, 

(7.1c) 

K c ijkl v k ,i Nj + G(ny = v v\ k a k a\Nj - qa^N, 

+ ng l on (0, T) x r 0 , 

(7-ld) 

«(V)e^(f!;» 3 ) 

a.e. in (0,T), 

(7. le) 

V = u 0 

on 11 x {t = 0 } 

' (7-If) 

where 

N(n) = -c mjkl [(n,m -V,j Smj)v\k 

],i in fig, 

(7.2a) 

G{n) = c m:,kl [(n,m-v,j-Smj)v\k)Ni on r 0 , 

(7.2b) 

and 

h l (t, •) = ~[c ijkl (u 0 + tw 1 + t —w 2 ) k ,i 

],j in fig, 

(7.3a) 

g l (t, ’) = [c Zlkl {U 0 + tW! + ^«l 2 )‘,l] 

Nj on r 0 . 

(7.3b) 

Solutions to 13.21) will be obtained as the limit as k — 

> 0 of solutions to 17.111. 

Suppose that v £ Wt is given. Let q = Id + / ' 

v and let al be the quantity 

U 

associated with 77 through 13.11). 
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We are concerned with the following time-dependent linear problem, whose fixed- 


point w = v provides a solution to (□): 

wl-u(a J l afw l ,k),j+( a iq),k= f°r] in (0,T) x , (7.4a) 

a k iW \ k = 0 in (0,T) x fl f 0 , (7.4b) 

w\ - 4c ijkl w k ,i],j +N(rj) i =f + nh l in (0 ,T) x Q s 0 , (7.4c) 

k c ijkl w k ,i Nj + G( V y = v w\ k a k ajNj 

- qajNj + Kg 1 on (0, T) x T 0 , (7.4d) 

w(t, •) G Hq (fi; R 3 ) a.e. in (0, T), (7.4e) 

w = mq on H x {f = 0}. (7.4f) 


Remark 5. The two forcing functions fUia\ ) are introduced for compatibility con¬ 
ditions at time t = 0, allowing the solution of m to satisfy (wt(0),wtt{0)) G 
I7Q(f2;]R 3 ) 2 and even to satisfy the same initial conditions as solutions to \d.2i) 
would. 

In the following, for the sake of notational convenience, we will denote by 

3 

N(uq, (wi)i =1 ) a generic smooth function depending only on Ei ||w3_i||_fU(ng;R3) + 

*=0 

||^ 3 -i||#;(nAR 3 ) ] (with the convention that uo = wo), by -lV((<7j)i=o) a generic 

2 

smooth function depending only on E 1192—i||ffi ( Q/;R 3 ) and by M(f,ng,Kh) a 

2=0 

generic smooth function depending only on ||/||v*(T) + ll/lkj»(T)+« [ ll«o||if 4 (ng ; R3) + 
lki||u 4 (ng;R3) + ||ai2 1| h 4 (Qq ;R 3 ) ]• Then, let w 3 G L 2 (R;R 3 ) be defined by 


w 3 = v [{a]a k w l 0) + (/o? ? )t t (0) - [ajq,j]tt( 0) in O 7 , (7.5a) 

= ftt(° ) + cmjkl [K V,m -V,] Smj)v,k],i] tt ( 0) in > ( 7 -5b) 

where the time derivatives are computed with any 77(0,2) = x, w = r/t(0) = uo, 
0 "MO) = w n , (n = 1 ,2 ), d?q(0 ) = q n , (n = 0,1 , 2 ). 

Let us now define 

b K (f>) = K(c vkl w 2 k ,i ,f,j )L= ( ng ; R), (7.6a) 

c K (<f) = K(c ljkl w^,i , 4>\j ) L 2 (n . ;R) , (7.6b) 

d K (4>) = n(c ljkl u 0 k ,i, )L2(ng;R) ■ (7.6c) 


By proceding as in [5|, we can establish the existence of a fixed-point to the 
system ED- This follows the lines of J5J by first approximating by a penalty 
scheme the divergence-free constraint in the fluid in our Lagrangian setting, and 
by performing a regularity analysis of the solution of fra allowing the use of 
the Tychonoff fixed-point theorem. Given the estimates obtained in [5], no new 
difficulty arise, since the parabolic artificial viscosity in the solid controls the forcing 
coming from the quasilinear part on a short time which is a priori shrinking to zero, 
and for this reason the proof is omitted here. 

This leads us to the following 

Lemma 7.1. There exists T K > 0 depending a priori on k and on a given expression 
of the type Nq(uo, (wij) 3 _i) + Nq (((7i)f =0 ) + Mo(f,K,g,K,h), so that there exists a 
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unique solution (w K ,q K ) £ Zt k of the regularized problem m- Moreover, w K £ 

v 4 (t k ). 


In the next section we will study the limit of these solutions of the smoothed 
problems as k —> 0 (this being problematic since the solutions to these regularized 
problems are a priori defined on a time interval shrinking to zero as k —> 0). 

Moreover, the following variational equations (for n = 0,1,2) are satisfied for 
any test function 0 £ L 2 (0, T K ; Hq (fi; R 3 )): 


J iPt 0)i 2 (f2;R 3 ) dt + V J (d t (OfeOfeWKir ), l 2 (Qq-,M . 3 ) 

+ K J 0 ( cijkld ™ WKk,l,< t ,i,j ' )L2 ( n °’ K l dt ~ J 0 ( d t( a k qK )’ 0 k ’l)L 2 (n';S) dt 

+ [ {c t 3 kl d?[{r),i-r},j - 6 i:j )ri,i ],<j>, k )y ( o-R3) dt 

Jo 

= J o ( d t F ^) L 2 {n f. R 3) + (dr/^L^fi-R 3 ) + dti^] b M) +a t n Wc K (</>) + d?[\]d K {4>) dt 

(7.7) 


together with the initial conditions w K ( 0 ) = uq, (w K )t( 0 ) = w i, (w K )tt( 0 ) = W 2 and 
g«(0) = qo, (g K )t(0) = qi, (q K )tt(0) = (? 2 - Moreover for the third time differentiated 
problem in time, we also have that a.e. in (0,T K ), 

2 II ( w K,)ttt (t)\\ l 2 (Q;R 3 ) v J (i a k a k w Kir)ttt, {w K )t.tt,s ).L 2 (fij( ;R 3 ) 

+ K f (c ljkl (w K ,)ttt k ,l (t),(w K )ttt\j 

Jo 

/ / \ < bi)tt [^{. a i)tt{ w K)tij F ^( a i)t( w K)tt>j F i a i)ttt w K >j}t 

Jo Jn 2 

+ / (<7 n)tt(t) [3(aj)tt(w K ) t y +3(a|)t (w K ) tt) j -\-{aP)ttt w K ij](0 
Jn f 0 


n [ 3 (a J i )tt(qK.)t(w Kj ) ttt ,j + 3 (a J i )t(q K ,)tt('w Kj ) ttt ,j +(Oi)ttt 9 «( ? ^«)ttt,j ] 

!q 

+ [ [ c tjkl [(r),i-r), j -5ij)ri,i]ttf {w K )ttuk 
Jo J ng 

< -ZV(u 0 , (Wi)Ll) + •^ r ((?i)i=o) + / [(^t«, ( w K.)ttt) L 2( n f. R 3) + (fttt, {w K )ttt)L 2 ( ng;R 3 )], 


(7.8) 


where, recall that C does not depend on the artificial viscosity k. The following re¬ 
sult will be fundamental to our proof that the time interval of existence of solutions 
to 17.111 is in fact ^-independent. 

Lemma 7.2. The mapping 7 : t —> || (w K , qK.)\\z t is continuous on [0, T K ], 

Proof. The continuity with respect to t of the terms of the type L 2 (0, t; H s ) is ob¬ 
vious, and since w K £ V 4 (T K ) (due to our artificial viscosity), so is the continuity of 
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En=0 


— n (Oq;M 3 )) * The only terms that remain are ||<9fw K ||L-(o,i;L 2 (n ; R 3 )) 


^3 

jn =0 I 

and ||5 t 2 g K || icxJ ( 0:t;i2 ( n / ;R ^. 

In order to treat them, we will invoke the fact that due to our artificial viscosity 
in the solid, we in fact have dfw K £ L 2 (0, T K ; L 2 (fl; R 3 )), which provides dfw K £ 
C°([0, T K \; L 2 (tt; R 3 )). For the second time derivative of the pressure, we notice that 
from the variational form, true almost everywhere on [0, T K \ for any 4> £ Hq (fl; R 3 ), 

(dt < / , )l 2 ( 0 ; R 3 ) T v{pt ( a k a k w Kir )i )L2 (q/.]j3) T k(c ^ ' d t W K ,; , (f) ,j )i 2 (Qg;R) 

— (9t( a k9 k) — i 4 1 >1 ) L 2 (Hq ;R) + ( c J [(Vii 'Vij <A)fc )i 2 (f2g;R 3 ) 

— (dtF,4>) L 2(fj^;R3) — {&t f i 0)L 2 (ng;R 3 ) — ^k(^) = ( a fc9«tti j)l 2 (Q^;M) > 

and the Lagrange multiplier Lemma 13 of |)] associated to the continuity results pre¬ 
viously established, we have the continuity of t —> II(o t-i 2 (n / -R)) on [OTk]- 

We now explain briefly why such a control on the fourth time derivative of 
w K holds, and is possible only with the addition of the artificial viscosity in the 
solid. In particular, this norm cannot be controlled as k —> 0, which is not crucial 
for our purposes in any case. In order to understand the idea, we return to the 
level of the setting of the fixed-point argument, where we assume that v in an 
appropriate convex set of Vf(T) x V^(T) is given, and search for a solution w of 
171 by a Galerkin approximation on a penalized problem (for the pressure), in 
a way similar to 0. The penalization parameter e > 0 is given, and we denote 

2 ' £71 J . 

— q n - -a{(w™) 1 ,j, where w™ is solution of the Galerkin approximation 


Qe = 

n =0 

at rank n, and where a\ is computed from 77 associated to the given v. Our interest 
will be with the first problem that appears in our methodology in ( 5 |; namely, the 
highest order time-differentiated problem is multiplied by d£w™ (which is permitted 
since it belongs to the appropriate finite dimensional space), and then integrate from 
0 to t. We obtain 

J ll 5 i 4 <lll 2 (n;R 3 ) + v j ( S (( fl M<.r),3t«) ) ») I2(fl / ;8 3 ) 


1 , L 2 (£l 


+ [^(c^5 3 (<) fc , I ,a 3 (<) i y ) i 2 ( o 5;R) ];-^ (a 3 (aUD, #«)*,. 

- [ (c ljkl d?[{r), i ■ri,j -5ij)r],i ], d 3 «U ) L 2 (ng ; R 3 ) 

= f[ [$F-d}w?+[ dU-dfw% 

Jo Jn s 0 J ng 

leading us for a time small enough depending on the artificial viscosity k (but not 
on n and e) to an inequality of the type, 


ll^t <llL 2 (n;R3) + sup[||5 t <|| H 1 (n /. R3) + n\\d t <lljri ( ng;R 3 ) + e ll 5 i <£\\ L 2 (n / ;B 

[0,t] 

< C e [N(uo, K) 3 = i) + JV(( 9i )?= 0) + N(f)], 
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where C e depends a priori on e. By proceding in a way inspired by our methodology 
in Section 9 of [5J, we can then prove that we have control, independently of e, on 
the first three norms. Taking the limit first as —> oo and then as e —> 0, indeed 
provides us with dfw K G L 2 (0, T K ; ; R 3 )) as announced. □ 

We note that this latter regularity property in the solid is only possible with the 
artificial viscosity k > 0. 


8. Estimate for the solutions of (17.41) independently of k 

In this section, we will denote ( w K ,q K ) = (w,q) and denote the corresponding 
quantities aj by aj. In what follows, <5 > 0 is a given positive number to be made 
precise later when it will be chosen to be sufficiently small. 


8.1. Energy estimate for w ttt independently of k. We are now going to use the 
regularity result (w,q) £ Zt k in the energy inequality 17.81) (which was established 
independently of the artificial viscosity), this time by interpolating and using the 
energy properties of the nonlinear elasticity operator, in order to get an estimate 
independent of the artificial viscosity. 


Step 1. Let I\ = 



C (Vii'VTj ^ij )R ! U -I 'IT tttik ■ 


An integration by parts in 


time shows that 


h 



'IT ttik 



'IT ttik (^) 5 


and thus with the properties of the Bochner integral in fL 2 (S4g;R), 


V,i ‘V,j (t) - Stj = / [fj,i ■w,j +w,i ], 
Jo 


we deduce 


\h \ <Ct ||lTtt|||«.(o,t;ffi(ng;R3))ll ? T||L”(0,t;/I 3 (Og;R3))||i7||L~(0,t;ff3(ng;R3)) 


<Ct llwl 


Wt- 

t 


Step 2. Let I 2 = 3 / / c ljkl (rj,i ■fj,j -6 l j) t w t ,i ■w tt t,k■ Similarly, 

Jo Jus. 


h = -3 


C [(l?;z 'Vij fiij'jt'Wttil 'ITtfifc ‘V-j 'IT ttik ] 


0 JUS 


+ 3 [f c ljkl (fj,i ■fj,j -8ij)tw t ,i ■w tt ,k (*)]o■ 
J US 


( 8 . 1 ) 
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By the same type of argument as for the previous step, we then get 

|/ 2 | < Ct ||w t t||ioo( 0)t;fl -i(n. ;R 3))||«)||L«'(0,t ; H3(n» ; R3))||?7||£,oo( 0)t;J j3(n. ;R 3)) 

+ Ct ||wtt||L«>(o,t;ff 1 (f2g;K 3 ))ll^lli“(o,t;ff 3 (n5;R 3 ))ll^' t ll-C'° 0 (o,t;ff 1 (f2g;K 3 )) 

+ Ct || vj t t\\L°°(o,t-,m(n‘-,R3)) IIw t |||oo(o it ; ff 2( 0 « ;R 3))||??||L~(o,t;ff3(ng;R3)) 

+ C ||c uW [(uou +Uo,i)w!,i + / -fi,j )t]||Loo(0 )t; z,2 ( ng;] K 3)) 

Jo 


x H^)fcMllL“(0,i;L 2 (ng;R 3 )) 


+ N(u 0 ,(wi)i = i)i 


and thus, 


\h\ < ^IHIwy +Ct\\w\\w t +CsN{u 0 ,{wi)i =1 ) . 


Step 3. Let / 3 = 3 / / c ljkl (fj,i -fj,j - Sij) tt w,i ■w tt t,k■ By an integration by parts 

Jo Jo* 

in time, 

I 3 — 3 / / C J [(^7,-i ‘T]-j fiij'jtt'UJtrf "'Mttik ‘Vij tl ’^ttik ] 

Jo Jtl‘ 


+ [ [ c ljk \v,i -f),j -8ij)ttw,i ■w tt ,k (*)]q- 


Similarly as before, we get 

|/ 3 | < Ct ||lhtt||L~(0,T;ffi(ng;R3))l|w|ll,~(0 > t;iL3(ng;R 3 ))ll^llL' !o (0,t;ff 1 (f28;R 3 )) 

+ Ct ||rt)tt|| L~(0,T;JL 1 (ng;R 3 ))||wt|||oo( 0it;ff 2( n s ;R 3))||??||L~(0,t;i/ 3 (ng;R 3 )) 

+ Ct ||w t t|||oo(o i T;Hi(ng;R 3 ))ll« : 'llL“(0,t;ff 3 (f2g;R 3 ))|l»?l|L“(0,t;ff 3 (ng;R 3 )) 

+ Ct ||lT«||L~(0,T;Hi(ng;R 3 ))||wt||L~(0,t;H 1 (^o; R3 ))ll'^llA 0 °(0,t;-f/ 3 (ng;R 3 )) 

+ C \\c ljkl ((fj, i -fj, j ) t t(0)w 1 ,i + f {{fj n ‘V,j —8ij)ttW,l )t)l|L~(0,t;L 2 (ng;R 3 )) 


x || U>tt,k ||L~(0,t;L 2 (OS;R 3 )) 


+ N(u 0 ,(wi)i =1 ), 


and therefore 


\h\ < $\\w\\w t + Ct \\w\\w t +C s N{u 0 ,(w l )i =1 ). 


Step 4. Let h = / c ljkl (fj,i ■fj,j -5ij) t ttf),i ■ Wttt,k • By symmetry of c, we notice 

Jo J% 

that 

-^4 — 77 [ [ ^ {Vn "Vij)ttt {Vil ’V>k)tttt 

^ Jo ./fig 

n c ljkl [4(fj,i ■w t t,k ) + 3 (rj,i -ri,j)ttt{wt,i ■w t ,k )], 

!» 
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and thus, 

V*-\ [ clJkl (v,i 'Vij)ttt(Vil ■V,k)ttt(t )I 
4 J ng 

< Ct ||Vw tt |||« )(0jt;L 2( n g ;R 9))||Vlij|| L oo( 0it;ff 2(ng;K9))||Vr?|| i cx ; ( 0it;/i -2( n g ;R 9)) 

+ Ct ||Vui4i oo (0,t;L 2 (ng;R 9 ))||Vwti|U°°(0,t;L 2 (ng;R 9 ))||Vw|||oo( 0jt;ff 2 (n g. R 9)) 

+ Ct ||Vii)t||| o =( 0i t ; ffi(ng ;R 9))||Vii)tt||L“( 0 ,t ; z, 2 (ng ;R 9))||Vj?||L“>( 0 ,t;fl-2(f2« ; R9)) 

+ Ct ||ViS|| L - (0 , t;ff 2 ( n S; R9 )) ||V«; t ||ic O(0)t;ff i (nS;R 9 ) ) + N(u 0 , (wi)^ =1 ) 

< N(u 0 , (wi)f =1 ) + Ct ||w||tv (8-4) 

Now, Let I = - j c ljkl (fj,i -fj:k)ttt(t). By expanding the integrand 

4 J ng 

with respect to the time derivatives and using the relation in iL 3 (fl§; R 3 ): fj(t, •) = 
Id + / w(t', •) dt' and estimates similar as in the previous steps, we find that 

Jo 

1 1- [ <J 3kl Wttn™ k tC (*)| <CsN(u 0 ,(wi)? =1 ) +8\\w\\%f t +Ct\\ui\\^ Vt . (8.5) 

Step 5. By using llOl and |H3 we find that 

| h~ [ c ijkl w J tt ,iW^ t ,i(t)\ <C 6 N{u 0 ,(wi)^ =1 )+6\\w\\^ t +Ct\\w\\$ Vt . (8.6) 

dng 

Step 6. By proceding in a way similar to 0, except that we replace the con¬ 
stants C(M) appearing there by appropriate powers of ||(iZ\ q)\\z t , we find that the 
integrals set in the fluid domain are bounded by 

S\\(.w, q )||| t + Cs{N((q z ) 2 i=0 ) + N(u 0 , (w»)- =1 ) + M(f, Kg, re/i) + C ||(w, q) |||j. 

Step 7. Thus, from (17.811 . and Steps 1-6, we then get on [0,T K ]: 

sup || u>t££ II x,a (n ;R 3) + j [ 11 djttt \\~ffl (q/ ;R 3 ) + K ll^t«||ffl(ng;R3)] + sup ||«)tt||fl-l(n»;R3) 

< Cs[N{u 0 ,(wi) 3 i=1 ) + M{f, Kg,Kh) + N{{qi)? =0 )+ti\\(w,q)\\% t ]+C6\\{w,q)\\% t . 

(8.7) 

Step 8. The estimate of qtt in L°°(L 2 ), independently of re, will require some adjust¬ 
ments with respect to the methodology of JSJ. To this end, we notice that we can 
apply a Lagrange multiplier Lemma similar to Lemma 13 of |^, but corresponding 
to the case a\ = 5j, to the variational form true on [0, T„]: for all <fi £ 7Lg(fl;R 3 ), 

(wttti JJ) L 2 (Q;R 3 ) + )tt , 4>,s )z/ 2 (njJ;R3) + k(c 2 ' W tt ,i , (f> ,j )i 2 (ng;R) 

+ (c ^ [{fjri 'Vtj ~ ]tt, J’lk )i 2 (ng;R3) ~ ((<^fc ~" ^kl)4tt, 4> iz).L 2 (fi£;R) 

— {JflkOjtt ~ 4* >1 )L 2 (n^;R) ~ i^tt, L 2 (Q. S 0 ;R 3 ) — (/it) < / > )L 2 (ng;R 3 ) — b K (<j>) 

= (qtt, div 0)i 2 (Q/. R p 
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which provides for any t £ [0,T K ]: 

II Qtt || L 2 (Hq ;R) — C[ll^tti|li 2 (n;R 3 ) + II )tt || L 2 (£2^;R 3 ) + K ll ^tt || H 1 (fig;R) 

+ WKVn -v,j -M^mIIl^r 3 ) + II {a l k q)tt - dl<7ti|U 2 (n ; R) 

+ ||a — Id|| II 94*11 _L 2 (fijf ; r) d" N(uq, {wi) j =1 ) + M(f , Kg, k h)]. 

By using <E3 for the first four terms of the right-hand side of this inequality and 
remembering that the L°°(0, t; L 2 (ff q ; R)) norm of qtt is part of the norm Z t for the 
next two terms of this inequality, we get 

sup||<Ztt||^ 2(n / ;R) < Cs[N(uo,(wi)? = 1 )+M(f,Kg,Kh) + lV((gi) 2 =0 )] 

[ 0 ,£] 

+ C s ti\\(w,q)\\% t +C6\\(w,q)\\ 2 z t - (8-8) 


8.2. Estimate on wtt and w t • From the previous estimates, and the arguments 
that we will see hereafter for the case of ui, we have 

H^ 4t llL 2 (0,t;7/ 2 (n^;R 3 )) II9tt IIL 2 (o,t ; ff 1 (Q f 0 ;R)) + INt|lz,~(0,t;if 2 (ng;R 3 )) 

< C s [N{u 0 , (u>i)i=i) + M(f, Kg , Kh) + N((qi)i =0 )] 

+ C S ti\\(w,q)\\% t +C6\\(w,q)\\ 2 z t - (8-9) 

Similarly, we infer from 18.911 that 

H^' t llL 2 (0,t;ff 3 (n^;R 3 )) H^llL 2 (0,t;/f 2 (n^;R)) 11 ^ 11 (0,i ;TT 3 (f2g ;M 3 )) 

< C s [N(uo, (Wi)f =1 ) + M(f, Kg, Kh) + 3V((gi) 2 =0 )] 

+ C 5 ti\\{w,q)\\ 6 Zt +CS\\(w,q)\\ 2 Zt . (8.10) 


8.3. Estimate on w. 

We will denote R(j_ = {x £ R 3 | X 3 > 0}, R 3 = {x £ R 3 | X 3 < 0}, and B_(0,r) = 
B{ 0, r) (~l R 3 . We denote by d/ an H 4 diffeomorphism from B( 0,1) into a neighbor¬ 
hood V of a point xo £ r 0 such that d>(.B(0,1) HR 3 ) = V nf2g, d/(£?(0,1) HR 3 ) = 
V H fig, d>(i?(0,1) n R 2 x {0}) = V (1 To, with det Vd* = 1. We consider a cut-off 
function £ compactly supported in 13(0,1), and equal to 1 in B{ 0, |). 

With the use of test functions <p p = — [p p * (£ 2 w o d/)], aiQ , lQ , 2Q , 2a3Q , 3 oil/ -1 (which 
is in L 2 (0, T k ; H) (fl; R 3 ))) in 17.711 for n = 0, and by denoting W = wo$, Q = god/, 
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E = fj o \]/, we get after integrating by parts appropriately and letting p —> oo, 

1 f* 

2 IICW,a 1 a 2 a 3 WIIl 2 (R 3 ;R 3 ) + 

] (Wt [C 2 W] 5a 1 a 2 a3 C ^ 0 : 10 : 20:3 )i 2 (R 3 ;R 3 ) 

~b ^ ]iaia 2 0!3 > [C W] jsaict 2 ct3 )i 2 (R 3 ;R 3 ) 


/O JR 3 . 


[Q fe il 7aiQ'2Q'3 [C 2 W? JJQ'IQ'2 0'3 

+ « r([C^ W (W,i-^-) *,iWa 2 a 3 . [C 2 Wq,fca ia 2 a 3 )l 2 (R5_;R 3 ) 

JO 

"b / ^ (E-i 'E.j Utji ■tEbj ) ^,i],aia2a3 ? [C BJ],fc aia2a3 )l 2 (R^;R 3 ) 

Jo 

< C N(uq, (wi) i= i) + lim / [ (F, 4>p) L2iQf.^3) + (/) ^p)L 2 (ng;R 3 )] 

p—t-oo Jq 'O’ 


f* t' 2 

J [ “b t C K (<^>p) + d K ((j> p ) ], 


lim 

p^°° Jo 


( 8 . 11 ) 


where C*** = c mn °Pg^gfgf,g l p £ R 3 (B(0,1);R), 5 = [V^]" 1 £ ff 3 (R(0,1);R 9 ), 

M = $(*)£■ 


Remark 6. Note that this limit process as p —> oo for the nonlinear elastic energy 
is possible because df"w £ L 2 (0, T K ; 7J 4_n (flg; K 3 )) (n = 0,1,) due to our artificial 
viscosity in the solid. Whereas we could also use difference quotients, it appears that 
the product rules are less cumbersome with the use of horizontal derivatives instead, 
which is permitted since we already know at this stage the regularity of w and q. 
Also, the limits on the right-hand side of is. in) do not present any difficulty, given 
the regularity of the forcing functions and three integrations by parts with respect 
to horizontal variables. 


Remark 7. Since £ is compactly supported in B(f), 1), the integrals set on R 3 , R 3 , 

R(j_ do not depend on the extension that we chose for W, E or Q, and simply 
represent a more convenient way to write these integrals. 

Step 1. Let L\ = K f ( [C 4 HI,; ],a 1 a 2 a 3 j [C W],kcna 2 ot 3 )l 2 (R 3 ;R 3 )- By 

Jo 

using the H 3 regularity of the coefficients C y , 

Ll = K f {C 9 W,i ai0l2a3 ’\E f,j \I/,J , () lR,feQ 1 a 2 «3 )l 2 (R 3 ;R 3 ) 

Jo 

+ K [ ( %i], aiaaaa A?,i, C 2 w, kaia2a3 ) L 2 (R 3 

JO 

+ « / ([C ijkl W,i -'*,,] 5 O 1 O 2 O 3 5 [C ^ // iifc]?0i0203 C 2 w, ka lft 2 a 3 )l 2 (R^;R 3 ) 

Jo 

>Ck (C ijkl W, ia lft 2 a 3 ‘1^0 ; C bL,fc QlQ , 2 Q , 3 )_L2( R 3_ ;R 3) Ck I II D)|| H 3 (Dq-R 3 ) 

Jo Jo 

(8.12) 


— Ck / ||VL,aia 2 a 3 II B 3 (B_ (0, i);R 3 ) 11 ZD 11 v 


;R 3 ) 
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Step 2. Let 


-^2 f ^ ^n ->j ) ] , 0 : 10 : 20:3 f [C -^t] ikai 0 - 20:3 )l 2 (E^_ ;R 3 ) • 

Jo 

£3 denoting the set of permutations of {1, 2, 3}, 

-^2 / {C ^ i^E ,2 -E,j Bila. 10 : 20:3 5 C E'tikatia20t3 )l 2 (K^_;R 3 ) 

Jo 

“I - 2 / (C C ^ (^nQift2a3 ) E,l 1 Et,kot.\a.2<y.3 )l 2 (R^_;K 3 ) 

Jo 

+ [ (C 2 [[C ijkl (E,i -E,j )], aia2as -2C iiW (£7,te 1 a 2 a 3 -E,j )]E,l , Et,ka 1 a 2 a 3 )_L 2 (R 3 ;R : 

JO 

E I ([C (C* 4 )jaia2C«3 E,l]t, E, kaiCt2 a 3 )l 2 (R 2 ;R3) 

Jo 

~ [(C (C* j )iaia2a3 E,i , -E' ) fcaia 2 a3 )l 2 (R^_;R 3 )] q 

^ ' f ([C [C 1 4 {E,i 'E,j — )] ) a < 7 (i) ■®'>o: < 7 ( 2 )a CT ( 3 )i ]iai t Et, a20 c 3 k )_L 2 (R^;R 3 ) 

<tGS 3 ^° 

^ ] / ([C [C 1 4 (-E'ji "E,j — ’F,* )])a<7(i)C« <T (2) ■®')a (r (3)/ ]iai i Et, a2 a 3 k ).L 2 (r5_;R 3 ) 

creS 3 J ° 

d" I {[C ^ {E,i "E,j \I/,j ) E,i ] , ai a 2 a 3 j [C ®i])aia 2 a 3 t — C Et,aia 2 a 3 k )_L 2 (R^;! 

Jo 

From the regularity of w and the if 4 regularity of \F, we then infer 

^2 = / (C 4 (E,i‘E,j 'Fjj ) E, aia 2 a 3 i, () Et, aia2Ce3 k )l 2 (r^_;r 3 ) 

Jo 

+ 2 / (C ^ E, aia2a3 i "E,j () E,i , Et, ai a 2 a 3 k )l 2 (R5_;R 3 ) 

Jo 

+ / (Ciail 2 a 3 E,i'E,j ( 2 E,i, E t , aiCi2Cl3k )i2( R 3_ ;R 3) + ^2) 

Jo 

with 

ILjI < <5||wi||u/ t + Ct ||w||^ t + CsN(uq, {wi)i =1 )- (8.13) 

By integrating by parts in time, we deduce 

1 r* 

L 2 = J (C 4 (E,i -E,j tF^ •\F,j )( E , ai a 2 oi 3 l , () E, ai a 2 a 3 k )i 2 (R^;R 3 ) 

+ (C 4 (E,i ‘E,j — d/,j ,j ) E, ai a 2a3 i , £ E, ai0l20 , 3k )l 2 (r 3 _;R 3 )] q 

2 / (C 4 E) aia20C3 i -E,j E t ,i , E , aia2a3k )^ 2 (r-^;r 3 ) 

Jo 

+ [(C* ^ E, aia2CX3 i ‘E,j () E,l , E, aiQ , 2a3k )^2fR3^ ;R3j] Q 
— / (Eia\a 2 a 3 (.Eii'E,j C )i; E, aiol2a3k )^ 2 (R 3 ;R 3 ) 

Jo 

+ [(C)ai k a 2 oc 3 E,i-E,j ( 2 E,i , E, ai0l20l3 k )l 2 (R 3 _ ;R 3 )] 0 + ^2) 
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which implies in turn 

\l J 2 ~ (C ^ C ^ 101020:31 (0 ' ^,2 1 P ia 1 a 2 0/ 3 k (i))i2(R3 ;R3) | 

< ^ll^llVt + CsN(uo, (« 2 *)i =1 ) + Ct ||w||v^ t - (8-14) 

With efc (A; = 1,2,3) denoting the canonical vectors of R 3 , let 

-P(A) = || 2c 4 ( fjymnil 'Cj') 6 fc — [c ^ " (?),i ■ fjj (5j_j ) Vik ] jlmn 11 L 2 (f!g-R 3 ) (^) ’ 

where to and n are arbitrarily fixed in {1, 2, 3}. We then have 

P(t) < Pi(t) + P 2 (t), 

with 

Pi (0 = || 2c4 {f),Umn'&ij ) Gfc 2[c " ) (fjimni ‘Vij )Vik ],2 11 £ 2 .r3) (A), 

P^it) = || [c J (fj,i-fjjj — 6ij)fj,k ],lmn +2[c J {fjimni ‘Vij )Vik]il || L 2 (f2g ;R 3 ) W ' 

We first notice that 


Pl[P) i^|| 2c4 ‘(j]rilmn 'Vij *)Vik ] 

+ || 2c J ’ [7),j mra ‘(f), j ] V,k ),l ||^2(fJg;R3)(A)- 


L 2 (f2g;R3)(*) 


Next, by writing 77 (f) = Id+ / ii) and [rj,mni ‘(v,j] Vik ),i ( t ) = [r), m m ] V ,k ),i (0)+ 

J 0 

/ [[»?,mm ‘{Vij ] Vik ),i ]t respectively in iJ 3 (fig; R 3 ) and L 2 (n s 0 ; R 3 ), we obtain 


Pi(t)<C[ |H| ff3(f2g;R3)] SU P[|| ? )llff4(ng;R3) + II 2? 11 // 4 (£2q ;R 3 ) ] 
J 0 [ 0 , 4 ] 

+ N(U 0 , (Wj)f = i) + C / ||w|| ff 3(ng;R3) sup ||7)||^-4( n s R 3) 
JQ [0,4] 

< N(u 0 ,(wi)? =1 )+Ct ||w||^ t . 


Next, we see that 

P‘2 (A) £ ||c^ [{Vii'Vij &ij)Vikmn P(Vii ‘Vij ),m Vikn ‘Vij )m Vikm ] ,2 || £2(qs. R 3) (A) 
+ ||c - [(?), m j ‘Vinj +?7,m ‘Vijm )fjik ],2 || i 2 (fJg;R3) (A)j 
and by the same type of arguments as for Pi(t), 

P 2 (t) < N(u 0 ,(wi)? =1 ) + Ct \\u>\\w t i 

implying 

P(t) < N(u 0 , (wi) f =1 ) + Ct ||w||^ t . (8.15) 

Now, from the definition of a solution to the smoothed problem ED, 

|| KC 4 (Vtiilmn ‘^b') ^2c + [c 4 ' ( Vii'Vj ^ij ^Vik ] imnl 11 ^2 , R 3) (A) 

= ||(^4 ^ ^ Kh), mn (t) ||d 2 (Oq ; R3), 

which implies with ED 

I mm ) + L{Vmm )||i 2 (ng;R 3 ) (t) 5; || (ult ~ f ~ Kh), n m (t) || L 2 (Og;R 3 ) 

+ N(uq, (rci)i =1 ) + Ct ||u , ||w t - 
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Since this inequality also holds for any t' € (0 ,t), Lemma KTT1 provides 

\\L(fj,nm )||L“(0,t;L 2 (ng;R3)) < C\\w t ,mn || L” (0,t;L 2 (fig;R3)) + C M (/, Kg, kK) 

+ N{u 0 ,(wi) 3 i=1 ) + Ct 

which with the estimate on w t from the previous subsection leads to 

II L(r},nm )IU-(0,t;L 2 (ng;R3)) <C s [N(u 0 , (wi) f =1 ) + M(f, Kg, Kh ) + N((qi)? =0 )] 

+ C s ti\\(w,q)\\% t +C6\\(w,q)\\ 2 z t - ( 8 - 16 ) 

Step 3. From the estimates on L 1 — L 2 , and similar estimates that we could get 
in the fluid as in [5], but this time by replacing C(M) by appropriate powers of 
||(w), q)\\z t , we then deduce that for all t £ [0,T], 

l 2 /'* 2 

^ (Oil L 2 (R3 ; r3) + V J (£ b k W, aia2a3r , W, Q , lCt2Ct3S )_L2(R3_;R3) 

< Cs[N(u o, (Wj)f =1 ) +M(/, Kg,Kh) + N((qi)? =0 )\ 

+ Cgt 4 || {w, q)\\% t + C5||(u;, q )||| t . 

By the trace theorem, we then get 

/ || W||^3.5 (S;R 3) <Cs[N(u 0 , (Wi) 3 = i) + M(f, Kg, Kh) + N(( qi )? =0 )} 

Jo 

+ CstH(w,q)\\ 6 Z t +CS\\(w,q)\\b 

where S = {{xi,x 2 ,xs) £ K 3 | |rci| < |rc 2 1 < \, X 3 = 0}. By a finite covering 
argument, we then get 

/ IMIiL3-5 (ro;K3) < C s [N(u o, (wO f =1 ) + M(f, K g, Kh) + N(( qi )?= 0 )] 

Jo 

+ Cst* ||(w, g)||| t + C(5||(w;, g)||| t . (8-17) 

From the estimate on wt and the trace estimate Em we infer in a way 

similar to by elliptic regularity arguments that 

H^L 2 (0,t;ff 4 (n£;R3)) + I ^11 L 2 (0,t;Lf 3 (n^;R)) 

< Cs[N(u 0 , K)?=i) + M(f, k 9 , Kh) + N(( qi )t 0 )| 

+ C 5 t*\\(w,q)\\% t + CS\\{w,q)\\ 2 Zt . (8.18) 

Similarly, from Em and the trace estimate (l8~T71) . elliptic regularity provides 
ll^ll!°°(o,t;ff 4 (ng;R3)) < Cs[N(uo, (wii)i=i) + M(f , Kg, Kh) + N((qi)? =0 )] 

+ Cst 4 \\(w,q)\\ 6 Zt +CS ||(w,g)||| t . (8.19) 

9. Time of existence independent of k 
From E3, EH), EH, EH3). and (l8~T9l) . we then have for any t £ [0, T K ], 
\\{w,q)\\ 2 Zt < C s [N 0 (u 0 , (wi) 3 ^) + M o(f, ng, Kh)) + N 0 ((qi)i =0 )\ 

+ C s ti\\(w,q)\\% t + C Q S \\{w,q)\\ 2 Zt . 

The subscripts 0 in Co, Nq, Mq mean that we no longer consider generic constants 
from now on. 
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Now, let Jo > 0 be such that CqSq = For k > 0 small enough, for any 
t £ (0, T k ) we have 

\\(w,q)\\ 2 Zt < 4 C' < 5 o [fV 0 ( UO) ( U ; i ) 3 = i) + Mo(/) + 7 Vo(fe)?=o)]+ 2 C 5o ti \\(w,q)\\% T , 

(9.1) 

where M 0 (f) = M 0 (f, 0,0). For conciseness, we will denote C\ = 2 Cg 0 and Ni = 

4C 5o [1Vo(wo, (^*)i=i) + M o (/) + JVo((®)?=o)]- 

Now for t £ (0, T k ) fixed, let a t (x) = x 3 -- H-—, so that 

C\t~i C\tz 

a t{\\{w,q)tz t ) > o. 

Now let t\ = [ 2 7 c^a ^ ] 4 > w hich does not depend on k, and let T = mm(T K ,ti). 
From now on, we assume that t £ (0,T). We then have a t {(2>C\t^)~^) < 0 which 
implies that a t has three real roots z i, z 2 > Z 3 , with Zi < — {3C\t*)~v < Z 2 < 

1 1 7 V"i 

( 2>C\tz )~2 < Z 3 . From the product Z 1 .z 2 .z 3 =-- and at(3N\) < 0, we infer that 

C±t 4 

0 < z 2 < 3iVi < Z 3 . From 19.111 and the continuity of t —> ||(to, ?)||z t (established 
in Lemma Q we then infer since ||(w), q )\\% 0 < -W < Z 3 that we have 

Vf £ (0, T], \\(w, q)\\ 2 Zt < Z 2 < 3Ni . (9.2) 

This implies that fj(T) £ tf 4 (f^; R 3 )rW 4 (fig; R 3 ), u)(T) £ i^(ft; K 3 )rW 3 (flg; R 3 )n 
H 3 (flg;R 3 ), w t {T) £ 7L 0 1 (L!;M 3 )nld 2 (L!g;R 3 )nld 2 (L!g;K 3 ), ih tt (f) £ H^fl-R 3 ), 
wttt(T) £ L 2 (fl;K 3 ), g(T) £ tf 2 (l^;R), q t (f) £ g«(T) £ L 2 (^;K), 

with a bound that depends only on the right-hand side of ED- The compatibility 
conditions for the smoothed problem o at T are also satisfied by definition of a 
solution, which means that we do not have any new term of the type of b K , c K or d K 
associated to w(T) to add to the already existing forcing terms coming from t = 0. 

We can thus build a solution of the smoothed problem 17.111 defined on [T, T+5T], 
ST depending solely on the right-hand side of iOl . that we will still denote {w,q). 
It is then readily seen that ( w,q ) £ Zf , ST and is a solution of the approximated 
problem o on [0,T + ST\. If T = ti, we have our solution defined on the k 
independent time interval [0, fi], with the n independent estimate 19.211 . Otherwise, 
if T < f 1 , we can also assume that T + ST < t \, which implies, in the same fashion 
as we got ED, 

Vt £ [ 0 , T + ST ], \\(w,q)\\% t <N 1 + C^ \\(w, q)\\% t - (9-3) 

This implies in turn that fj(T+6T), w(T+5T), wt(T+8T),wtt(T+ST), wttt(T+ST ), 
q(T + ST), q t (T + ST), qttif + ST) are in the same spaces as their respective coun¬ 
terparts at time T, with the same bound as well, since we could from E3 repeat 
the same argument leading to 19.21) . this time on [0,T + ST]. The compatibility 
conditions at T + ST being also automatically satisfied, we can thus build a solu¬ 
tion of the approximated problem 17.111 defined on [T + ST,T + 2 ST], the time of 
existence being the same as starting from T from the similarity of the bound that 
we obtain on fj(T + ST), d?w(T + 5T){n = 0,1, 2, 3 ),d?q(T + 5T)(n = 0,1, 2) and 
their respective counterparts at time T. We will still denote this solution ( w,q ). 
It is then readily seen that ( w,q) £ Zf +2 sr an( l is a solution of the approximated 
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problem on [0,T + 2 ST], We then have in the same fashion as we got (19.111 . 

Vte [0,f + 2ST], \\(w,q)\\% 1 <JV 1 +C' 1 ti ||(«i,g)||| t . 

By induction, we then see that we get a solution ( w,q ) defined on [0,ti],satisfying 
the estimate 

Vt€ [ 0,h], \\{w,q)\\z t < 3iVi = 12C 5o [iVo( W o,(^)i=i) + Mo(/) + A r o(fe)?=o)] ; 

(9.4) 

establishing the independence of the time of existence respectively to k, since t\ 
does not depend on k. In the following we will note T = t\. 


10. Existence for (Id. 21) 

Proof. We can here choose to take k = i, and let n —> oo. By the bound 
independent of n on [0, T], we then have the existence of a weakly convergent 
subsequence of (u>, q) in the reflexive Hilbert space Yr, to a limit that we call (v, q ), 
which also belongs to Zt and satisfies the estimate 

\\(v,q)\\z T < 3 Ni = 12C 5 o [N 0 (u 0 , (Wi)f =1 ) + M 0 (f) + N(( qi )? = 0 )]. 

The usual compactness theorems ensure at this stage that (v, q) is a solution of (Id. 21) 
on [0, T], The smoothness of our solution ensures that the solids do not collide with 
each other (if there is more than one) or the boundary (for an eventually smaller 
time), which establishes the existence part of Theorem. 15,11 □ 


11. Uniqueness for (Id.21) 

Proof. Since we cannot use a contractive mapping scheme for our problem, we have 
to establish uniqueness separately. Let then (v, q) denote another solution of (Id.21) 
in Zt- Then, taking v — v as a test function in the variational formulation of 
the difference between the systems (Id.21) associated to each solution provides for 
te [0,T]: 

1 T f 4 

2 II( u — WIIl^h-rs) + v {a r k a s k v, r -a r k a s k v, r , v, a -v, a ) L 2 (n /. R 3 ) 

+ / (c ljkl [{-6 ij)fj,k], V,l-V,l) L 2 ( ng;R3) 

Jo 


Jo ( a i q - a i q ’ vi ’ j - vl ’j )j 


->/.»! = J (/ °V- for),v-v) 


L 2 (Oq ;R 3 ) 


( 11 . 1 ) 


13 nj 

For the viscous term in the fluid, we write 

a r k a s k v, r —a k a k v, r = a r k a s k (v, r -v, r ) + (a k a s k - d r k a s k )v, r , 

which with the L°°(0, T; H 3 (fl q ; M. 3 )) control of v and v provides us with an estimate 
of the type (where C denotes once again a generic constant) 


J { a k a k v ir a k a k v iri v is ^is )/,2 . jj 3 ) — C j 


\v — v n2 


m( nj ; R3) 


-C 


( 11 . 2 ) 


0 Jo 


1 AT 1 (fi^;R 3 )' 
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Concerning the forcing term in the fluid, we first notice that if we still denote 
E(Cl)(f) as /, 

f(t, fj{t , x)) - f(t, rj(t , x)) = f f,i (t, (r) + t'(fj - r]))(t, x))dt' ( ff(t , x) - rf(t , x)), 

Jo 

leading us to 


II fit, 7? (t, •)) - f(t, v(t, -))ll L1.5(^;E3) 

< Cllr?^,-)-r?(C-)llL6(n^;R3)E / f x))dxdt']° 5 , 

with t, x) = rj(t, x)+t'(fj(t, x) — rj(t , &)). We have </>(f', t, •) € R 3 )(~lC 1 (On 

rg; M. 3 ). Moreover <j>(t',t,dCl) = dCl. We then have by invariance by homotopy of 
the Brouwer degree (for the parameter t) 

Vz G 0, deg(^(f', t, •), 0, 2 ) = deg(0(t', 0, •), Cl, 2 ) = deg(Id, Cl, z) = 1, 

which together with the regularity of <j>(t', t, •) establishes that <j>(t', t, -)(f2) = Cl and 
that Cardl^ -1 ^', t, -)(x)} = 1 for almost all x G Cl. Thus, 


la 


f,i (t,<t>(t',t,x))dx = / f,i (t,y)\detV(t>(t',t,(t> 1 (t',t,y ))| 1 dy, 

f J <p{t' 


which with the L°°(0, T; H 4 (CIq; R 3 )) control of ?y and ry provides 

[ fn <£(*', t, x))dx < C [ /,? (f, 

J J Q, 

Consequently, 

||/(t, fj(t, •)) - f{t, r){t, -))|| L i.s in f. R3) < C\\ fj(t, •) - r)(t, -)|| H i (a / ;M 3)||/||i/i(a;M3), 
implying 

| [ (forj-f °fj,v-v) L2{n , m \ <C , ^||/||i 2 (o it;ff i(n^ 3 ))||t;-w|| 2 2 ( 0 )t;irl(niiR 3 )) . 
0 

(11.3) 


Concerning the elastic term, 

{c tjkl [(v,i -V,j ~ 5ij)v,k ~{v,i-fj,j Sij)fj,k], v,i -v,i )L2 (n g. R 3) = /| + h + / 3 , 

with 

h= (c zjkl {v,i-V,j ~V,k), v,i-v,i) L 2 ( ng ; R3) 

Jo 

= -^(c ljkl (v,i -V,j Sij){v,k ~fj,k ), V,i-f),i )L3 ( ng;R3)(t) 

(c l3kl (r],i ■r), j -6ij) t (v,k-fj,k ), V ,1 ~V,i )l 2 (ng ; R3) 

— — ^II^W — *7(^)ll/f 1 (ng;R3) ~ C I Wv — ^llifi(fJg;R3)! 

0 

where we have used the L°°(0,T; iJ 3 (flg;R 3 )) control of u and v for the inequality. 
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Next, for the same reasons, 


h = I (c l3kl (.rin-Vn) ■ V,j v,k, v,i-v,i) L 2 ( ng;i 


= / (c l3kl (v,i ~V,i ) • (v,j ~V,j ) V,k , V,i -V,i )L 2 (n g ;i 


+ / {c tjkl (v,i ~V,i ) • V,j V,k, V,l-V,l) L 2 {n s R3) 

Jo 

= / - fj,i ) • ( T},j -fj,j ) fj,k , V,l -V,l )L2 (n g ;R3) 

Jo 

+ ) • V,j fj,k , V,i -V,1 )L2(ng;R3)(i) 

- / (c ljkl (r),i —fj,i ) • fj,j v,k, i?,;-?7.;)L 2 (ng;R3). 

Jo 

We then write for the second term on the right-hand side of the last equality 

fj,i (C') = + / v,i, to get by Korn’s inequality 

Jo 

h > C[\\q(t) — fy(i)||/fi(n» ; R3) — \\v(t) ~ f7(^)lli 2 (ng;R 3 )] — C^sup ^ — ^llff 1 (ng;R 3 )- 
Similarly, 

h= (c ufei (??,j -f},j ) • fj,i fj, k , V,l -V,l )L2 (n g ;R 3) 

Jo 

> C[ \\rj(t) — ??(i)||lfi(n»;R3) — \\q(t) — ^7(^)IIi 2 (ng ; R 3 )] — C [ Wv ~ ''?ll7f 1 (ng;R 3 )- 


Thus 

C 
/ 0 


/ {c ljkl Kvri -V,j -5ij)r),k-(v,i -11,3 -6ij)r),k], v,i ~v,i )z, 2 (fi5 ; R 3 ) 

Jo 

'i C[ WviO ~ viO lllr 1 (ag;R 3 ) ~ ll ? 7(0 ' ? 7(0lll 2 (ng;R 3 )] ~ C f \\v ~ ? 7llffi 

0 


(fig;R 3 ) • 

(11.4) 


Concerning the pressure term, with a 3 q — a 3 q = (a 3 — al)q + a\{q — q) and the 
L°°(0, T; ff 2 (flg; R)) control of the pressure, we get 

- J o ( a U~ afav 1 * -v\j )i, 2 (n / ;R) 

> -Cfv^llg - 9ll L oo( 0it;L 2( n / ;R ))||L’ - ^Ilz,2( 0it;ff l(f2/ ;R 3)) + t\\v - ^lll,2( 0]t;H l(Q/ ;R 3))]- 

(11.5) 

In order to get the estimate of q — q in L 2 (flg; K.), we have to introduce the time 
differentiated problem. By taking vt — Vt in the variational formulation associated 
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to the difference between the time differentiated systems, we obtain 

1 C 

2 I ( v t ~ ^t)(0lli 2 (n;R 3 ) "b V J ([ a k a k v ir ~0>k^k^fr]t> bjs — Vis ]i)_L2(Q/ ;R 3) 

+ / (c Z:>kl [{V,i-V,j- S ij)V,k-{'n,i-'n,j-Sij)fj,k]t, [v,l ~V,l ]t)L*((l‘-S*) 

Jo 

- J ([o<9 - [v\j -v l ,j }t) L 2 {n f. R) = J ([/ ° V - f ° v]t,v t - vd^nf.pay 

( 11 . 6 ) 

For the fluid viscous term, we easily find with the L 2 (0, T; if 3 (flg; K 3 )) control of 
the first time derivative of the velocity that 

J {[ a k a k v ir ~0'k®k'V,r]t, Vt is ~ Vt is ) £2 (q/ ; r3 ) > C{\ — t) J || V% — Vt || H x ;K 3y 

°(11.7) 

Concerning the forcing term in the fluid, since (/ o q) t = ( ft. + v l f,i )(rj) (with a 
similar formula for c), we then deduce in a way similar to the steps leading to m 
that 

rt 


| j ([/ ° il - f ° v]t, [v - v]t) L 2 (n f. t 


< CVt [||/t|Ua(0,t;fl-l(n;R3)) + ||/||L2(0,t ; if2(n ; R3))] ||v t - ^*11^2(0 t ■ H 1 (Q f -R3)) ■ 

°(ll-8) 

For the elastic term, we can also essentially reproduce the arguments leading to 

(trot . leading us to 

/ ( c l]kl [(v,i -V,j Sij)v,k -fj,j Sij)fj,k ]t, [v,i -v,i ]t)L2 (n « ; R3 ) 

Jo 

> C[ ||v(t) - v(t) ||ffl(ng;R3) - II v(t) ~ H(i)|||2 ( ng;R3)] - Ct sup ||t> - t>||ffl ( fig;R3). 

(11.9) 

The pressure term will require more care since we want to avoid the introduction 
of q t — qt that the most direct method would lead to. To do so, we notice that 

/ ([ a i9 — ®i9]t) b ij —V ij ]t)i2(n^;R) = ^4 + I 5 + ^6) 


with 


h = Jo ^ a ^ tq ~ 

h = j ( a {(qt — qt )> b ij —v ij ]t)z,2(f2/.R)> 


^6 — / ([af a\]qt, [v ,j v ,j ]t) i 2 (Q/ ;R )- 

J 0 

For Z4, we have in a way similar to 111.511 . 

\h\ < C[Vt\\q - q\\ L oc ( 0 )t;I ,2 (a /. m ^|bt - v t " 


L 2 (0,t;,F/ 1 (fi£;R 3 )) 

+ t\\v t - ^l! i 2( 0 it;ff i (n / ;R 3))]- 
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For J 6 , the L 2 (0, T; ff 2 (flg; R)) control of q t provides us with 


For I 5 we have: 


l^el — ^ Ut llL 2 (0,t;ff 1 (n^;R3))' 


^5 — J (Qt Qt,0, J i V t ,j a i v tij )L 2 J (?* 9ti ( a i a?K,i) La(n /;R) 

= J o (Qt - qt, (a J i)tv\j ~(di) t v\j ) L , (n / ;R) - J (qt - qt, (a? - d?)vl,j ) £2(n / ;R) , 

where we have used the relations ajv 1 ,j = 0 = d\v l in fig for the first integral. By 
integrating by parts in time, 

h= [ (Q-Q,[(ai)tv\ j -(dl) t v\ j ]t) L 2 {n f. m + f (q - q, [(a? - a-Kb }t) L 2 (n f. R) 

J 0 ^ 0 

+ (q-Q, ( aj) t v x ,j -(dj) t v\j ) L 2 (n f. m (t) + (q-q, (a? - dj)vlj ) L 2 (n f. R) (t). 
With the L 2 ( 0 , T; Ff 3 (flg; R 3 )) control of v t we have 

I J (q-q, \( a i)t v ,j ~(d{)td ,j ]t)i 2 (n^ ; R) 1 + 1 J (q — q, [( a i)t — I 

< eVt II q ~ 9llLoo(g it;i 2( n / ;R ))||r’i - ^t|lL 2 ( 0 l t ; Hl(^ ; R 3 ))> 

I (q ~~ Q, ( a i)tt> ,j~(d{)td ,j )i,2(n^;R)(^)| 

— CVt \\q(t) — g(f)|| L 2 (n / ;R) ||^i — ^i|lia(g it;if i(f 1 / ;R 3 ))* 
The remaining terms are more delicate. We first have 

I j Q (q~ 9, (o< - «<) )i 2 (n';R)l + |(« - 9, (Oj ^ 0() «ty ) L 2 (fi;(;R) W| 

— ^ J o Il9 — 9llL2(Q/ ;R )||a — «llL4(n 0 / ;R 9)l|V%|| i 4 K; R 9 ) 

+ lk(i) - 9(i)|| L 2 (n^;R)ll a ( i ) a WllL 4 (n^;R 9 ) l|V« t (i)|| 

L 4 (Oq ;R 9 ) * 

( 11 . 10 ) 


The apparent problem here is that a —a is estimated in T 2 (flg;R 9 ) in terms of v — v 
in H 1 ( flg;R 3 ). Now, a bound of this quantity in L 4 (flg;R 9 ) will require a bound 
of v — v in 7J 2 (flg; R 3 ). In order to get such an estimate, we will bound v — v in 
H 2 (i !q;R 3 ) by lower order terms in v — v. To do so, let us first estimate the trace 
of v — v on r 0 by using the test function — [£ 2 (v — v) o \F|, aa oil / -1 in the difference 
between the variational problems satisfied by v and v. By proceeding as in Section 
cm we would then get an estimate of the type, where S > 0 is given: 

J IICV[(v -v)o T], a ||i2 (R 3_ ;R9) + ||CV[(J7 -fj) o T], a (t)||i2 (R 3_. R 9 } 

<C[Vt + 5\ [||v-w||^ 2(n / ;H 3 ) + ll9-9||^i (n / ;R) ] +Cs J Q ll u « -^lli 2 (n ; R3) 

; R 3)+C/ ll'f - H|||fl(n;R3) 

J 0 


C 


ll ? ? vWh 2 ^* 
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which by patching all the charts defining To leads to an estimate of v — v in 
L 2 ( 0 , f; 7? 1 - 5 (ro; R 3 )) yielding by elliptic regularity: 

ft 


W- v\\1„ f + \\q-q\\ 


'ff 2 (nj ;R 3 ) 
ft 




< C[\/f + (5] [||v - v||^ 2 (n / ;R3) + Ik - 9||ffi (n / ;R )] + Cs J o Ik* “ ^tIIz, 2 ( o ;R 3) 

+ C/ 11*7 “ *?lllt 2 (ng;R 3 ) + C / 

Jo Jo 


Ik Olif^fi;! 


Thus, with a choice of S > 0 small enough, we have for t small enough and the use 
of Gronwall’s inequality, 

IkW _ *ik)ll/r 2 (ng;R3) + [ Ilk - ^Hff2 (n / ;R 3) + Ik - £||^i (fJ / ;R) ] 

J 0 

< C 1 / Ikt-^tll| 2 m ;R 3 )+G / |k-u||^i (n . R 3 ). 

Jo Jo 

By using this estimate in mint . we then get for a time small enough 
I f k - T (O* - «;) ««« )z, 2 (Qq ; r)I + |(9-9, (ffl< - 0<) ^td )i 2 (n*;R)k)| 

<cVt[J lkt-^t||l 2 (n;R3) + J Ik _ ^llffi(n ; R 3 ) + Ik - 9l|^oo(o^ ;L 2( n / ;R ))]- 

By putting together the estimates on Li, Is and Ig, we have 

| f ([ a iQ — ®i9]t, k ij ij ]t)i2(fj/ ;R ) I 


< cVt [ f Ikt - «tll£ 2 (n;R3) + / Ik - + Ik — 9llj 

Jo Jo 


”'I.L'»(0,t;L 2 (fi i *;R))]' 

( 11 . 11 ) 

Now, by considering the difference between the two variational forms satisfied repec- 
tively by (v,q) and (v,q), and writing the difference between the pressure terms 
as 


[ (aU~aU)4> l ,j= [ al{q~q)(t>\j+ [ (aj - a])qf ,j , 

Jn f 0 J% Jsv 0 

the Lagrange multiplier Lemma 13 of 0 provides for all t € [0,T]: 

Ikk) — 9(011 L 2 (f2^;R) — II ( v t ~ ^t)(OII L 2 (Q;R 3 ) + II ( v ~ v) (0 II H 1 (£2j(;R 3 ) 

+ II (*7 — ^?)(^)II Af 1 (f2g;M 3 ) + Veil'll — H|k2( 0jt;i ji(f2/ ;R 3))]- 

( 11 . 12 ) 

By putting together the estimates irra - irmi . we then obtain for t u > 0 small 
enough an inequality of the type: 

Ikt _ ^t|ll~(0,t u ;L 2 (n ; R 3 )) + / Ikt _ ^t||^i(Q/.R3) + Ik — ^lll°°(0,t„;/t 1 (f2g;R 3 )) ^ k 

which shows that (v,q) = (v,q) on [0,t u ]. Let 

T u = sup {t e [0, T]| (v,q) = (v,q) on [0 ,t]}. 
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If T u < T, we can repeat the same procedure with T u replacing 0, which would lead 
to uniqueness for [T u , T u + St) as well. Thus, we have T u = T, which concludes the 
proof of the theorem. □ 

12. Optimal regularity on the initial data 

We first remind some extensions and regularization results on domains: 

Lemma 12.1. Let f Y be a domain of class H A . Then, there exists a linear and 
continuous operator E(Ll') from R 3 ) into 7L m (R 3 ; R 3 ) (for each 0 < m < 4) 

such that E(Cl')(u) = u in f Y. Also, if the H 4 norms of a family of domains stay 
bounded, the norms of the corresponding linear operators also stay bounded. 

Lemma 12.2. Since fig is of class H 4 , let if m £ if 4 (i?_(0,1);R 3 ) (m = 1 
be a collection of charts defining a neighborhood of its boundary. We note 

N 

H\\ H 4 — J2 H^HI H 4 (B- (0,1);M 3 ) • 
m =1 

Then, there exists a sequence of domains (fig",) of class C°° , so that fig C fig’” , and 
which are defined with a collection of charts if m ’ n £ R 4 (13_(0,1); R 3 ) (m = 1, ...,N) 
so that Y ) m =l HV’ m - h 4 (B-( o,i);R 3 ) —^► 0 as n —> oo. We then denote the 

complementary of fig" in fi by fig'” and Tg = <9fig”. We also assume n large 
enough so that the different connected components of fig” ( if there is more than one 
solid) do not intersect each other or the boundary of ft,. We denote a n = ||fl s,n ||#6. 

We now state the optimal regularity assumptions needed in our analysis, and 
explain the adjustements required to the previous proofs. 

Theorem 12.1. With the same assumptions as in Theorem fO except for the 
following concerning the regularity of the initial data: 

u 0 £ lL 6 (flg ; R 3 ) D H 3 (flg; R 3 ) D -ffg (fl; R 3 ) D L 2 div j, (12.1a) 

/s(0) £ lL 2 (flg; R 3 ) (~l if 3 ' 5 (r 0 ; R 3 ) , (/ s ) t (0) £ if 1 (fig; R 3 ), (/ s ) tt (0) £ L 2 (fl; R 3 ), 

(12.1b) 

the conclusion of Theorem EH holds. 

Remark 8. We have chosen here to take different forcings for the fluid, that we 
still note f with the same assumptions as in Theorem \5.1\ and the solid, in order to 
stress out that the higher order regularity required indeed comes from the hyperbolic 
scaling of the Navier-Stokes equations. The somewhat not so natural condition 
f s ( 0) £ Lf 3 ’ 5 (Tg;R 3 ) is made in order to get w\ £ iJ 4 (flg;R 3 ) associated to the 
condition W 3 £ L 2 (flg;R 3 ). 

Proof. The idea is to first regularize the domains, initial data and modify the 
forcings in an appropriate way, and then pass to the limit. 

Given 0 < p £ V(B(0, 1)) with f B ^ Q ^ p = 1, we define as usual p n (x) = n 3 p(nx). 

f 

We first notice that uq, wi, qo and q\ still have the same regularity in fig as in 
Theorem 15.11 We first define in fig’", itg = uq and w™ = w\, q(f = qo, qf = qi 


30 


D. COUTAND AND S. SHKOLLER 


(which is permitted since fig’™ C fig. We next define w% in fig’™, 

—vAw% + VgJ — Pn* E{ fig)(—+ Vq2) in fig'", (12.2a) 

div ii >2 = -[{ai)tt(0)uo,j+2(al) t {0)wl,j] in fijj’™,. (12.2b) 

w% = 0 on <9fi, (12.2c) 

V dW^- ( & Nn = + -p n *E(n f 0 )(q 2 )N n on rj, (12.2d) 

where TV" denotes the unit normal exterior to fig’™, and finally Wg £ L 2 ( fig’™;R 3 ) 

by 

u>3 = W{ai4u, k ),j-{a 3 i q, j )f =1 + F]tt(0) in fi 

where the time derivatives on the right-hand side are computed with the usual rules 
from u(0) = uft, d p u{ 0) = w” {p = 1, 2), d p q( 0) = q* (p = 0,1, 2). 

We next define in the solid Ug by 

L 2 Uq = L 2 [p n -k E(Qq)(uq)] in fig ™, (12.3a) 

K = «) f on rj, (12.3b) 

LK)+p„*E(fiS)((/ s ) t (0)) = K) j on rj, (12.3c) 

where the right-hand sides of the previous boundary conditions come from the fluid 
regularization previously carried. Note also that 

Lu%£H 4 { fi*’”;R 3 ), (12.4) 

(with an estimate that may blow up as n —> oo) since 

L(Lu5) = L 2 [p n kE{n s 0 )(u 0 )] in fi*’™, 

L(uS) = -p„*E(Q s 0 )((f s ) t (0)) + K)^ on T™. 

We can then define fg in fig” by 

L 2 /o = L 2 [p n *E(n°)(f s (0))] in fi*’™, 

/o™ = (w?)* on T™, 

C mjki (/g", m -Id,, +fS,j -Id, m )N£ = -2 c mjki «, m -Id,,- +<„■ -Id, m K,i 

+ i/ [u*,fe afaf] tt (0)lV™ - [gra^] tt (0)lV™ on Tg , 

with the same conventions as for the previous system for the time derivatives eval¬ 
uated from fig’™ and c m ^ kl (ug , m -Id,,, +Ug ,j -Id, m )ug N™ evaluated from fig’™. 
We then define in fig™, 

< = /o", 

u >2 = [[c mjfci (r?, m -??,i -Smj)v,k},l] t (0) + p„*£(fig)((/ s ) t ( 0 )), 

= LK)+p B *£(fia((/ 8 ) t (0)), 
w 3 = [[c r ™^ w (r?, TO - 7 ?, i -( 5 m j)??,fc],i] tt ( 0 ) + p n *£(fig)((/ s ) tt (0)), 

where the time derivatives on the right-hand side are evaluated with u(0) = Ug, 
Vt{ 0) = u>™. We also define the regularized forcing in the solid 

/"(*) = p n kE(n s 0 )(f s (t) - fs(0 )) + /o”in fi*’™. 
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We then have Uq, w w% in JZ£(fi; R 3 )rW 4 (fl£’ n ; R 3 )n£T 4 (flg’ n ;R 3 ) anddivitg' = 0 
in flg’ re , W 3 £ L 2 (fl;R 3 ), with 

II^UK) - u o|| J? 4 (n / ;Ill 3 ) + ||i«")W) - Mo||ff3(n-R3) -> 0 as n -► 00 , 

(12.5a) 

II^UK) - W lllff4 (n / ;R 3) + ||£(nS’")«) - Wl|| fl -2(ng;R3) 

+ 11^2 - w 2 ||ffi(n;R3) + ||B(n^ n )K) - w 2 || ff2(n / ;R3) -^Oasn->oo, (12.5b) 

\\ u 0 l|ff 6 (n“’";R3) < Pn, ||^r||if4(n“'" ; R3) < Pn, 11^2 l|ff 4 (n“'”;R3) < Pn, (12.5c) 

11^3 - w 3 || i 2 (n . R 3 ) -> oo as n -> oo, (12.5d) 

where p n is a given polynomial expression of a n and n. We briefly explain how 
those constants appear. For instance, for the first estimate of 112.5cll . we have by 
elliptic regularity on <m that \\uq ||ff6(f2®.»;R3) is bounded by a sum of terms, 
one of which being P(||fig "||#e) 11(^2 )^llff 4 (n / ’"- R 3)’ E being a polynomial which 
does not depend on n. Next, still by elliptic regularity on (t TF3 . we have that 
||(w 2 )* II fl-4 (n /,n ;R 3) is bounded by a sum of terms such as ||p„*P(fl6)(Aw2)|| ff2 ( n /,*v R 3) 
This particular term, by the properties of the convolution, is in turn bounded by 
n 3 ||.E(f!Q’")(w2)||ffi(R3 ;R 3). This shows that a term of the type P(a„)n 3 ||w 2 || H i( n / ;R 3 ) 
appears in the sum of all terms bounding \\uq ||ir 6 (n s '”- ;H 3)- Since the other terms in 
the sum can be dealt with similarly, this explains our estimate (ES3- 

For the pressures, we have 

\\m f o’ )(Qo ) _ 9o||/ i -3(o/. R ) + ||£ , (^o )(?") ~ ( Zilli/3(n/. R ) 

+ \\ E (^o n )(Q 2 ) - 92 || J ji (n /. R) -> 0 as n -> oo. 

( 12 . 6 ) 


Since the initial data Uq and forcings /"(0), /"(0), / t "(0) are smooth enough to 
ensure the regularity properties dim we then deduce that we have similarly as 
for theorem ED the existence of a solution w n of a system similar to o with 

f 

/, uq, fig, fig being replaced by their counterparts with an exponent n , and b K , c K , 

d K being replaced by b n , c n , d n (with the choice k = —— --) given by 

n{pn + 1) 


bnify = 


c n {4>) = 


1 


n(p n + 1 ) 

1 


f c ijkl on n k 


n k li \ 

2 i! ,j lZ, 2 (n*'";R)) 


,ijkl n k 
W 1 >Z ) 


ft ij ) L 2 (f2p’ n ;R) + ( w 2 ~ 


n{p n + l) v v "° 

+ 0)-^"]? =1 + [(afaf)u,fe] t (0) JVj\ 


- (c t3kl [(T},i-T),j -^j)%]i(0)JVr, ^)L 2 (r";R3), 


w 2 , < / > )l 2 (q/.™ ;R ) 

0)l 2 (T£;R3) 


d n(<t>) = n(/ J + y , <t>\j )i>(nS--;R) 

+ ( _ [( a i9)(0)N J ”] 3 =1 + [(a^af)u,fe](0)lV", (/Oz^r^R 3 ), 

where the time derivatives are computed with a velocity satisfying iz(0) = u q, 
ut(0) = u>i and a pressure such that q( 0) = qfi , q t { 0) = qf. Note that by construc¬ 
tion, the solutions w n to these problems in fl satisfy w n {0 ) = Uq, w nt (0 ) = w", 
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Wntt{ 0) = Wnttti 0) = w 3 ■ Next, we proceed to energy estimates similar to Sec- 

tion|Hl The bounds obtained are similar, except that this time the terms associated 
to b n , c n and d n tend to zero as n —> oo. This is clear from the convergence results 
C23J, imm for the integral terms associated to the fluid. The terms associated to 
the solid asymptotically tend to zero by properties of the convolution. For instance, 
with the notations of Sectional for <p p = — \p p * (C 2 Wn ° ’F)],aiaia 2 a 2 a 3 a 3 ofl> , we 
get after change of variables, an integration by parts in time, three integrations by 
parts in space: 


t'- 


dtt'\ 


1 Jo 2n(/3„ + 1) ( 2 ,l ’ ^ p,j )i2 ( n o 

c 

- n (/3 , + 1 ) l|w 2 llff 4 (n“'”;R 3 )ll 7 7 ™||L~( 0 ,t;ff 4 (n^";R 3 ))) 
and thus with our estimate 112.5cll . we have 


where Z™ denotes the same type of space as Z t with fig and fig being replaced 
by their counterparts with an exponent n. This type of estimates thus shows 
that this term does not change the energy inequalities in Section |H1 We can thus 
reproduce the arguments of Section |£l establishing that (w n ,q n ) can be defined 
over a time T independently of n, and that its norm in ZJJ depends solely on 
N(u g, {w™)i = 1 )+N((qi)? = 0 )+M(f n , 0,0) and thus, thanks to the estimates (112 5all . 
(I12.5bll , 112.5dll . solely on N(u 0 , (wi)? =1 ) + N((qi)? =0 ) + M(/, 0,0). We can then 
consider the sequence (F?(fi)(w n ), £(fig’ n )(g n )) which is bounded in a space similar 
as Zt, but defined on R 3 and extract (with respect to n) a weakly convergent 
sequence in a space modified from Yt by replacing the condition u £ Hq (fi; R 3 ) by 
u £ H 1 ( fi;R 3 ) . By the classical compactness results, we next see that the weak 
limit ( v,q ) £ Zt and is a solution of (Id.2D with / as forcing and u(0) = uq. This 
solution is also unique in Zt- D 


2 n(p n + 1 ) 


(d 


,ijkl n k 
w 2 il 


) L 2 (Hq’™ ;R) dt' 


< — 1 | 

n 


13. The case of incompressible elasticity 

In this Section, we explain how to treat the supplementary difficulties appearing 
when the incompressibility constraint is added in the solid. This leads to the same 
system as 13.211 . with the addition of the condition det V ?7 = 1 a. e. in fig and the 
addition of [(a^q),k ]f =1 on the left-hand side of (I3.2dl) and the addition of —qa^Nj 
(the trace of q being from the solid phase in this new term) on the left-hand side 
of I3.2el) . We now state our result and explain how to overcome the additional 
difficulties related to this constraint. 

We first update our functional frameworks. While Xt and Wt do not change, 
Yt and Zt become respectively 

Y t = {(v,q) £ X t xL 2 (0, T ; L 2 (fi; R))| d?q £ L 2 (0,T; H 3 ~ n (n f 0 -R)), 
d?q £ L 2 { 0, T; H 3 ~ n ( fig; M))(n = 0,1, 2)}, 

Z T = {(v,q) £ W T x L 2 (0,T;L 2 (fi;R))| d?q £ L 2 (0,T; 7? 3 -"(fi^;R)), 

dtQ € L 2 (0,T; U 3_rl (fig;R))(n = 0,1,2)| q tt £ L°°(0, T; i 2 (fi; R))} . 
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Remark 9. Whereas the pressure in the solid satisfies dfq £ L°° (0, T ; if 3_n (Rg; R)) 
(■ n = 0,1,2), it appears that the limit pressures q K are controlled uniformly in the 
norm of Zt and seemingly not in these norms. Note also that whereas the velocity 
field is smoother in the fluid phase for the solution of our next theorem, the pres¬ 
sure field is actually smoother in the solid phase. Whereas our artificial viscosity 
smoothes the velocity field in the solid, it also interestingly makes the pressure in 
the solid for the regularized system less smooth than the one associated to the solu¬ 
tion of the constrained problem, which is source of difficulties that we shall describe 
later. 

We now state our result: 

Theorem 13.1. With the same regularity assumptions as in Theorem II 2. 1\ and, 
assuming that the compatibility conditions associated to our new system at t = 0 
hold (for the sake of conciseness we do not state them here), the conclusion of 
Theorem fO holds for the case where the incompressibility constraint is added to 
the solid part. Furthermore, dfq £ L°°(0, T; fi 3 _n (Rg; R)) (n = 0,1,2). 

Proof. The extra regularity (with respect to the norm of Zt) on the pressure in 
the solid simply comes from the equation 

v t -c m3kl [(T), m -T], j -6i j )T], k ],i+a J i q, j = f in (0,T) x fig, 

which once the regularity for the solution w £ Wt is known provides immediately 
the result. We now explain how to obtain a solution in Zt- 

The beginning of the proof follows the same lines as for the compressible elastic¬ 
ity case. We first assume that the initial data satisfies the regularity assumptions 
of Theorem EH and define the same smoothed problem as o with the corre¬ 
sponding updates for the incompressibility constraint. We then define the same 
fixed point linear problem as EH where the condition affw l ,k= 0 in Rg is added 
(the a* being computed from the given v) and add a^q,k on the left-hand side of 
E3H and —qofNj (the traces being taken from Rg) on the left-hand side of itOdli . 

We then proceed as in [H] to construct a solution to this system by a penalty 
method (the penalty term being this time defined over R) and get the same type 
of regularity result. This provides us with a solution (w K , q K ), that we also denote 
by ( w,q ), for the incompressible version of 17.111 on a time T K shrinking to zero. 
As for the compressible case, (w K ,q K ) is in Zt k , and since our smoothed problem 
has a parabolic artificial viscosity, we also have for the velocity in the solid the 
regularity dfw £ L 2 (0, T K ; H i ~ n (fil g; R 3 )) (n = 0,1,2, 3) (with estimates that blow 
up as k —» 0). Thus, ( w K ,q K ) £ Zt k with 

Z t = {(«;,«) £ Z t | dfw £ L 2 (0,f;R 4 -"(Rg;R 3 ))(n = 0,1,2)}, 
endowed with the norm 

\\iw,q)\\\ = ||Kg)||| t \\ d ? w \\ 2 LHo,t-,H*-«(n‘-,R3)y 

n =0 

We next proceed as in Section 0 to get energy estimates, that will be carried this 
time for the n dependent norm of Z t , independently of n on [0,T K ], and for such 
a purpose it is important to keep the k 2 factor in the definition of the norm. We 
could extend the sum to n = 3, whereas it is not necessary. 

As before, the first set of estimates has to be carried on the highest order time 
derivative. Our energy inequality EH has the same form, except that the integrals 
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over Qq where q appears has to be taken this time on fh The part over S2 q is 
estimated as before. We now explain how to deal with the integrals set on fig for 
the pressure, which indeed needs some justifications given that the velocity in the 
solid is not controlled uniformly in n in a space as smooth as the velocity in the 
fluid, while the pressure is controlled in the same type of spaces in both phases. 


13.1. Estimates on uittt • Here t denotes any time in (0, T„). The most difficult 
integrals set in [0, t] x fig and associated to the incompressibility constraint in the 


solid are K i = 


o dng 


qttiaDtwlu,] and K 2 = 


o dn; 


q(al)tttwl tt ,j , the others 


being either less difficult or similar to estimate. 
Step 1. For K i, if we denote N s = —TV, we have 


\Ki\ = h 


/o 

r t 


( Qtt),j ( a^tw'ut + 


0 JTq 


Qtt(a 3 i)twl tt N? I 


< 


C[ / l|9tt||i? 1 (OS;R)||^||i? 3 (Og;R 3 )||w||j? 3 (nS; I R 3 )||^ttt||L 2 (OS;R 3 ) 

do 


+ 


m i 




I 3 (fig ;R 3 ) IMI 3 (Og ;R 3 ) 11 i 11 


H?( fig;R 3 ) J 


< CVt\\(w,q)\\% t +Ctz sup[ ||gtt||22 (ng . R) ||w t M || r 2 


[o,t] 


' L 2 (Qq ;B 


(W,l 




Zv 


(13.1) 

(13.2) 


where we have used the continuity of Wttt in the sense of traces along To to bound 
the L 2 (To; K 3 ) norm of Wttt by means of the K 3 ) norm. Note that we have 

also used the fact that the L°°(L 2 ) norm of qtt is in the definition of the norm of 
Z t . In order to get an estimate on this norm, we would proceed in a way similar as 
to get (EHJ) in Section El 

Step 2. Concerning K 2 , we have by integrating by parts in space 

K 2 = - [ f q, j (a J i )utwl tt + [ [ q(a > i ) t ttwl tt Nj , 

J o dog do dr 0 

since our artificial viscosity provides the regularity w t t € L 2 ( 0 , T K ; iJ 2 (flg; R 3 )) 
and Wttt G T 2 (0, T K ; R 3 )) (with estimates that may blow up as k —> 0). 

The difficulty here comes from the second integral. Whereas as for A'i we can 
estimate the trace of wttt on To from the fluid, we have to take the norm of S7w t t 
in Ff~ 0 - 5 (To;R 9 ), which is problematic given that the norm Z t contains only its 
L 2 (flg;R 9 ) norm. In order to circumvent this difficulty, we notice that the same 
formula holds if we replace w t tt by E(VlI){w{ u ) (the extension to R 3 of the velocity 
in the fluid). Since Wttt = u>{tt on To, we have Wtu = E(£Iq)(w{ u ) on To, which 
implies: 

K 2 = - [ f q,j (■aDutWttt + [ [ q^tttmi){Mu):) 

do dng do dng 

+ [ [ Qd(^i)tttE(n{,)(w{ tt y, 

do dng 
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and thus, 


K 2 < C 1 / l|9l|ff 3 (n;;R)ll^tt||ff 1 (^o; R3 )ll^ll' H ’ 3 ( n o; R3 )ll^' Mt ll i2 ( 0 ; R3 ) 

Jo 

Jo 

+ C Jo lift + J 0 9t||H2(n» ; R)||^ti||/fi(n«;R3)||»?||H3(n«;R3)||«; tt t|| JJ -i (n / ;1R 3) 

+ C J Iko + J 9t||_H-2(ng;R)l|Wt|lff 1 (ng;R 3 )||w|| fl -3( n g. R 3)||u; tM || fl - 1(n / ;R 3 ) 

< C \\{w,q)\\z t J lkl|iL3(n» ;R ) + C \\(w,q)\\% t \\qo\\HZ(n°-R) J ll^tttlli/^n^RS) 
+ eVt \\{w,q)\\ 2 z t \\qt\\L 2 ( 0 ,t-,H 2 (n°;WL)) J ||||jyl(j-j/;R3) 

< eVt [ \\{w,q)Wz t +N((qi)i= o) ]■ 


The most difficult integral set at time t on fig and containing q is 


I < 3 


qtt(aj) t wl t , 


for which we apparently just have an estimate of the type |Z 3 1 < C || (w, q)\\% (with¬ 
out any small parameter in front). We now explain how to treat this difficulty. 
Step 3. We first notice that 


K-i 


Qtuj {aDt.wlt 


'To 


qtt(a 3 i) t wl t N? 


If we could say that q tt is L°°(If 1 ) controlled, the L°°(L 2 ) control of w t t would give 
us a suitable bound for K 3 . Whereas we have seen in the statement of our theorem 
that qtt for the limit solution is indeed in L°°(H l ), we cannot seemingly get such a 
bound on the approximate pressures qtt- In order to get around this, we introduce 
similarly as in the previous step the extension to the solid domain of the velocity 
in the fluid. Since a similar integration by parts formula holds when we replace Wu 
b y mlKrtt), we deduce 


I < 3 


Qtt,j 


qtt(ai) t E(n f 0 )(w{ t y 


Qtt,j {dDt.Ein^iw^y. 


(13.3) 
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The easier term to estimate is K\ = [ q tt (ai) t E(n f 0 )(w{ t y tj , for which we have 

Jw 0 

for an arbitrary <5 > 0: 

\K%\ < C'||9tt||z,2(fJg ; R)||Id + J «i|| ff 3(fJg;K3)||Uo + J ■W" ; t|||f2( n 8 ;R 3)||fo||| r 3(Qg. R 3) 

x [ll«' 2 || if i (n / ;R 3 ) + Vt\\w tu \\ ^(o.t^Tn^jR 3 ))] 

< C , || 5 «IU 3 ( n- R) [l + i||(w,g)||z t ][Ar(u 0 ,(w; i )f = i) +t*\\{w,q)\\l t ] ||(w,g)||| t 

x [N(uo,(wi)^ =1 ) +t?\\(w,q)\\z t ] 

< C\\(w,q)\\ Zt [ N (uo,{wi)i=i) +ti\\{w,q)\\l t ]\\{w,q)\\l t 

x [N(uo,(wi)f =1 ) + ti\\(w,q)\\z t } 2 

< mw,q)\\% t + C 5 N(u 0 ,(w i % 1 ) +Cst*\\(w,q)\\% t . (13.4) 

For the hrst integral, the nonlinear elastodynamics equation in fig provides 

Vqtt = a~ l [-w m + nLwu + c lJkl [(v,i ■f ),j - 6 ij)fj,k]tt,l -2a t Vg t - attVq + ftt + nh], 
leading us for K\ = / q tt ,j {aDt.wlt to (since a -1 = Vrj in virtue of det Vrj = 1), 

K 3 = [ Y^V [-wttt + c zjkl [(rj,i ■fj,j -Sij)rj,k ]tt,i -2a t Vq t - a t t\7q + ftt + Kh]] J (al) t wl t 


k [Vfj [Lwtt]\ J (a 3 i ) t w t u . 
Jq.° 


(13.5) 


The integrals on the hrst line of this equality do not give any trouble and can be 
estimated in the same fashion. For instance, we have for 

K l = [ Nv[c ljkl [(v,i -V,j }] J (aj) t w l tt , 


\K 3 \< C\\ Id + / fo|| ff 3(Qj» ;R 3)||u’t||/i-2( n » ;R 3)||«;2 + / u)ttt||L 2 (ng ; R3) 

Jo Jo 

f t _ 1 3 

X I|lt0+ / 11 f[2 (Qs ; R 3) II fo || , R 3) 

J 0 

< ^\\(u’,Q)\\z t + Cst* W(w,q)\\ 7 Zt + C s N(u 0 ,(w i )^ =t )). (13.6) 

Now, the difficult term to handle is K 3 =n [X7rj[Lw t t]Y(a{)twl t . We hrst write 
the divergence form Lw^ t = cr,™ p ( wtt ), and integrate by parts: 

Kl = -K [ [Vfj,m [a mp (mt)]UYW)Mt 


-K [Vf][a mp {wtt)]l=iY[(ai)twl t \,rn 


+ K / m<r mp (*tt)]Ux} j (H)MtNn, 
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leading us to 

\ K i - K / [V? 7 [o- mp (w tt )]p =i P(d^) t «) t l t iV^| < CK\\{w,q)\\% t , (13.7) 

d r 0 

and thus by putting together (113.511 . 113.(ill and IH7I) . 

\[ Qtt,j (al) t w l tt ~ k f [Vfj[a mp (wtt)]l =1 Y{dl) t wl t N^\ 
dog J r 0 

< CK\\(w,q)\\ 5 Zt +C s ti\\(w,q)\\ 7 Zt + 5\\{w, q)\\ 2 Zt (13.8) 

+ Cs[N(u 0 , (m)i= i) + iV((gz)i=o) + M (/, «/»)]• (13.9) 


Now, the apparent problem comes from the term a mp (w t t) on To that should be 
taken in iJ _0 ' 5 (r 0 ; R), which is troublesome since the norm in Z t appropriate for 
our limit process only contains its L°°(0, f; R)) norm. In order to circumvent 

this, we notice that we also have, since E(ttl)(wtt) is at least as smooth as Wtt in 

fig, 

| / Qa,j (a\)tE(Q. f 0 ){wtt) 1 - k f [Vr)[cr mp (w tt )]p=i] J (Oi)t £; (llo)(^/ t ) I lV m | 
dog Jr 0 

< CK\\{w,q)f Zt +C s ti\\{w,q)\\ 7 Zt + S\\(w,q)\\% t 

+ Cs[N(u 0 ,(w i )i =1 ) + N((qi)i = 0 ) + M(f, Kg, /eft)], (13.10) 

leading us, since w = E(i lg)(w;l) on T 0 , to 

\( {al)twl t - [ q t t,j (al)t^(f2o)Kt) l | 
dog dog 

< C , «||(*,Q ; )||z t + C,5^||(w,g)||z t + <5||(w,g)||| t 

+ C s [N(u 0 , (wi)Li) + ^((9i)?=o) + M (/> (13.11) 

Thus, by using (troii . (trot) and (tom we have 

|/ 3 | <(C/c + C' (5 t*)||(u;,g)||^ +5||(w,g)||| t 

+ C 5 [iV(uo, (wi)f = i) + N((qi)i =0 ) + M(f , Kg, k/i)]. (13.12) 

Thus, we finally arrive to estimates analogous to <E3 and (lOli . with the right-hand 
side being of the same type as in 113.121) . 

13.2. Estimate on w t t and w t . With the same arguments as in the next subsec¬ 
tion, we have for n = 2,1: 

ll^t ^llL2(0 tt; fJ4-"(n£;R3)) + ||<9 t 9|| i 2( 0)t . J J3-n(nJ ; R3)) + ||^ t fj \ \ L aa (Q T . H 4-n (Q3 ;R3)) 

+ ^ II L 2 (0,£;7/ 4 - n (Og;R 3 )) + II ^11 L 2 (0,i;/7 3-n (Og;R 3 )) 

< C S [N(u 0 , (Wi)f =1 ) + M(f, Kg , nh) + N(( qi )? =0 )] 

+ (Ck + Cst*)\\(w,q)\\ 7 Zt + C8\\{w,q)\\ 2 Zt . (13.13) 

We now explain on the case of the highest space derivative how to obtain elliptic 
estimates independent of k , since the addition of the pressure term does not allow 
to use Lemma rb.il directly in the present case. 


13.3. Estimate on fj in fig. 
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13.3.1. Regularity of the trace of fj. . First, by proceeding as in Section |3 
and as for the case of the highest order time derivative, we get an estimate for the 
trace similar to (18171) . with a majorant of the same type as in 1113.121 . We explain 
hereafter how to handle the estimates related to the pressure in the solid in order 
to get this trace estimate since difficulties different than in the higher order time 
derivative problem appear in the higher order space derivative problem. 

Step 1. Let Qi= / [Q bj] 50 !iQ; 2 Q !3 [C 2 w% aia2a3j . 

Jo Jr 3 _ 

Then, 


Q i — Q 2 + Qo + Q 4, 


with 


Q2 I I Q bJ^ ai ot2a3 [C ] Tai<X2Ct 3 j > 

Jo Jr 3 _ 

Q 3 = / / [[ Q bJ ] ,0,10,2013 Qia ia 2 a 3 Q , ai a 2 a 3 ] [C W' ],aic 

Jo j Ri 

pt 


Qa= / Q, ai a 2 a 3 H[<?W'], ai 


0 J R 3 


For Q 2, we first notice that for 0 = fj o <F, if ( l i k is the sign of the permutation 
between {i,j, k} and {1, 2, 3} if i,j, k are distinct, and is set to zero otherwise, 


j T/T ri _ _ mni pqj [o m o n ] YAT i 

'i lOtiOi 2 CX3 VV lOL\OL 2 OL^j — fc I u iq \1aLiat2a3 vv Ij 

m e n w i 

ipa 1 a 2 &3 ’<? ijotia 2 a3 


2^ w l " ip ~ iq J50:10:20:3 ** ijatioc 2 0:3 

— f rnni f pqjn m f) n W i 

c c w ipot.1 0:20:3 w iq rr Jj 

1 1 \ ' e rnni pqj n m an rsr i 

W _ ^ c c ’pa„d) "i"” ”> 


crGE 3 


j 0 W 

’P“<r(l) 2)«<r(3) ’JOHCX 2 Ots 


1 1 V c mni f P < iba m a n uy* 

” r 9 Z_^ ’P“»(1)««(2) ’ga„(3) Kf !iaia 2 a3 

<xG£ 3 

_ pqj la m ni ]Q n 

2 C c r ^0:10:2^3 ’j a 10:20:3 

4. 1 W f mni pqj a m an Trri 

^ 2 ’P a &( 1) , <?«»(2)a 0 .( 3 ) rr ’jaia 2 a 3 

a G£ 3 

I 1 \ ' e mni pqj a m an Tjy * 

^ 9 ’P“<r(i)a<r( 2 ) ’9«ff(3) b“ia 2 a 3 > 


crGE 3 
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where we have use e mnl e pqq = e nm ''e 9W on the second equality and e mnl e pqq = 
e mm e jqp on the third one. Thus, 


Q2= QK, aia2a3 [[CW}\ aia2Ct3j -c 2 w; jaia2Ct3 } 


0 J R 3 _ 

l f f £ 

ae - S3 

1 


^ e mni e™[Q9™ ai 
<res 3 


f) n du W 1 

’p«r(i)0»(!) U, qa*( 3 ) hai KK U'a 2 C «3 

■pq 0 /- 2 n m ni !qq n l 

’paict2C>:3 ’jaiQ!2Q:3 ">q 

_J_ [ f (inni jpqj s-2 n m ni " 1 * 

o L / fc t ’ c, ’pai«2a3 (7 V“i«2a3 y °i u iq Jo’ 

^ JR 3 


showing 

|Q 2 | < CVi||(w,( 


|4 

II 


e mni e pqj^2Q 


K ai a aaa Q0,) (t) + N« qi )t o). 


’po;ia2C»:3 ’gai 0 : 20:3 ^ ">J 


In order to estimate the remaining term, we notice by integrating by parts twice 
for Q 5 = [ e mni ePqjc 2 e ^ ia2a3 e- qaia2a3 Q9 ,* that 


/-) _ mni pqj I m r\ n r\ r\ i 

V5 t t / S u ’qa 1 O 2 O 3 u ipctia2&3 ^ 


J x 3=0 

Since e pqq = —e qpq , we then infer 


hd 2 nd m 9 m 9 ” — i 9 m 9 n 

LVS ' £ v t7 ?j /)P 17 5 a 10:2 <23 ^’goi0:20:3 \S> )">q u "xx^ 0.2^3 u ’>pcx\<X2(X^ . 

C 2 Q 0 Z a 2 a 3 [^ 0.03 ( e 3 )p - 0.^203 M,]- 


2Q 5 = - e m ™e^ / 0 ™ lCt2a3 [(C 2 O0,j ), P fl,J aia3a8 -(C 2 Q0,} )„ 0,” ai a 2a3 ] 


+ e mn ’e MJ / C 2 Q 0,£ Q2 a, P 

J 2 : 3=0 


Jqoi O 2 O 3 ’’J 


e,)(e 3 ) p - 9,; aia2a3 0,5 teU 


Now, if we note 9? = E(Ql)(ijf) o \E>, we also have for 


/ mni pqj/-2nf m nn 

I c c S 17 5paio:2Q:3 17 iqcx\cx.2 03 Ko v l, i 
JR 3 


a similar formula. Since 0 ,™ lQ , 2Q , 3 = 0f 

2 Q 6 =-e mni e™ [ {0-9 f )£ iaaaa 

JR^ 

+ 2 f e mui ePqJ( .0 


” ia2a3 on {a ; 3 = 0 }, we then have 

kc 2 o 0 ,} ), p e, n qaia2a3 -(<; 2 Qo,i )„ ^ Q1Q2 « 3 ] 


1 n Of) 1 
’50:102^3 5 
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leading us to 


IQsl <C1|Id+ / w||^3(ng;R3)||f)|||f4 (n g ;R 3)||gO + / gt||j?2(f2» ; R)||Id+ / w|| ff 3(Og;R3) 


■C'||Id+ / ^ll^ 3 (fl /. R 3 ) ll^llH4 (f 2 3. R 3)lko + I gt||ff 2 (n g ;R) ||Id- 


^ 11 -f /3 (Oq ; R 3 ) 
r t 


+ C||Id + J , w> / H^r4(f2^ ;]R 3)||77||^r4 ( - f 2g ; j R 3) 11^0 + J 9t||/4 2 (ng ;R )||Id + J w|| fl -3(ng ; R3) 

< ^ll(^ I 9)lll t +CVt\\(w,q)\\% t +C s [N(u 0 ,(w i % 1 ) + N{(q i )^ 0 )] 

Step 2. We see by integrating by parts with respect to the direction a\ that we 
have 


\Q 3 \<CVi\\(w,< 

Step 3. Next, = Qe + Q 7 , where 


111 


Qg = 


/ 0 JR 3 _ 
rt 


Q-,a 1 a 2 ct 3 bjC W , 


0 JR 3 


Q 50:10:20:3 billowy 50:10:20:3^ C 50:10:20:3^ ]• 


Qi = 

We first have 

IQrl < C / [||g||_f/3( n g ;R )||ui|| ff 3( n g ;R 3)||^||| / 3 (n g ;R 3 ) ] < CVt\\(w,q)\\z t 
J 0 

For Qq the divergence condition b\W l ^ — 0 on Supp £ implies 

2 


^50:1020:3 C 2 [ Kw l 50:10:20:3,7 - Ww * ,j )iaiQ2a 3 ], 


Qe — 

Jo -tRi 
which in turn provides, 

I <96 | < C f [||Q||H3 ( ng;R)||w||_f/3 ( ng;R3)||^||^4 (n ,. R 3)] < C\ft \ \ (w, (?) \ | , 

J 0 

which concludes the estimates on the pressure terms in the solid, justifying why we 
obtain a trace estimate similar as with a majorant of the type of the right- 

hand side of (113.121) . Now, we turn our attention to the recovery of the regularity in 
the solid, which will need some justifications since we cannot directly apply Lemma 
lb.II 

13.3.2. Regularity in the incompressible solid. First, with the introduction of 
F = - ui t + c^ kl [(fj,i ■fj,j —6ij)fj,k],i -Lfj + f + nh in fig. 


and of f, the solution in fig of 


-f t + r = q , 


r( 0 ) = g 0 , 

we have for the nonlinear elastodynamics 


— Lw — Vrj Lfj + V[^r t + f] = Viy F in fig, 
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i.e., 


^[—Vi? Lfj + Vr]t — Vfj Lfj + Vr = V?) F — Lfj in fig. (13.14) 

We now apply Lemma IQ to this equation, leading us to 

sup || - Vi? L(fj) + < sup ||Vi) F - Lrj\\ H 2 <n‘■,«.*) 

[<M] [o,t] 2 

+ || — -t'(Id) + Vgo|| ff 2( n g ;ffi 3), 


and, with H = —L(fj) + X7r, to 


SUP ||i?||_f/2(ng;R3) 

[ 0 ,t] 

< sup[||V 77 F — -Vw Lr?||H 2 (n g ;R 3 ) + ||(V 77 — Id) L(f?)||ff 2 (n; ; R 3 )] +N((q i )1 =0 ). 

[o,t] 2 

(13.15) 


We then want to use elliptic regularity on the system: 

-Lfj + Vf = H, in fig, (13.16a) 

div 77 = (—a\ + 5ij)rj\j +3, in fig, (13.16b) 

fj = ? 7 |r 0 on F 0 , (13.16c) 


where the trace on Tg is estimated as we explained in the previous subsection. Now, 
for the divergence condition in fig, we notice that: 




= a i 
+ 


i «1*2 *3 V 
^ r~ n 


+(di - 5 i:i )fj l 


3 


frGS 3 


l a iji<r(l) V lji< t(2)*<t(3) ~^ a i V U*<r(3) 


For the apparently problematic first term on the right-hand side, we first notice 
that 


~J . fii . — _ F mni pqj(~ m~n\ ~i mni pqj ~ m ~ n i 
a ini 'I 13 ~ fc vhp Iiq Ml l hj ~ t t 'hpi 1 'hq '113 

— inm j qp ~ i ~n~m_ 9 %3 .. 

~ t 'hji 1 'hq 'lip — ^ a i'l 131 1 J 

which with the condition d J i fj' l ,j = 3, provides 0 = aj,^ f \ l , and thus 

®ijiii2«3 V 13 = —®inii2 V >ji 3 ~®iii 1*3 V iji 2 ~ iti V\ji2i3 • 


We then deduce that 


ll(«i l|fl-3 ( na; R )(t) < <5||(w,g)||| f +C s N(u 0 ,{wi)i=i) + Ct \\(w,q)\\z t - 

(13.17) 

Now, with 113.151) and (113.1711 . elliptic regularity on (113.161) provides for 

ll ? ?llz,“(0,t;ff 4 (ng;R 3 )) + ||f ~ t|| L”(0,t;ff3(ng;R)) 

a bound of the same type as the right-hand side of ifrrm with however the norms 
in Z t being replaced by the norms in Z t , due to the term k||ViZ> L? 7 ||# 2 (qs. R 3 ) 
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appearing on the right-hand side of firm that we bound by 

C«;||'u)(t)|| ff 3(Qg; R 3 )||I/(Id) + / iw|| ff 2( f2 g ;R 3) < C , K||li)(t)|| fl '3(Qg ;R 3)V / t||w||L 2 (0,t;ff 4 (ng;R3)) 

Jo 

< cVt\\{wA)\\%- 

We now turn our attention to the pressure, that we just need to control in 
L 2 (0,t;ff 3 (flg;R)). In order to do so, we notice from (lid.1411 that we have for 
K = f[-LtS + Vf t ]: 

K = ^ [V77 — Id] Lw + \ 7 fj Lfj — Vf + Vfj F — —Vui Lfj, 

which with the previous estimate on fj and f shows that we have a bound on 
||A'|| 2 2 (q t-H 2 (n 3 -R 3 )) °f the same type as the right-hand side of (113.121) . but where 
the norms in Z t are replaced by norms in Z t due to the estimate in L 2 (H 2 ) of 
k[V ?7 — Id] Lw. Now, elliptic regularity on the system: 

— Lku> + V nrt = 2 K, in fig, 

divKic = n(—al + Sij)u> l ,j , in fig, 

kw = kw, To on Tg, 
provides after integrating in time an estimate for 

K2 [||^lll 2 (0,t;ff 4 (ng;R3)) + ||ft - 1^1 ft 111,2(0,t;JJ3(ng ;R ))] 

with a bound similar as in firm still with the norms in Z t being replaced by 
norms in Z t . 

Thus, we obtain for ||g — / 0 „ g|| 2 2 ( 0 t -H 3 (n s -«)) the same type of estimate as 

well. Given our estimate on q t t, this also implies the same type of majoration for 

ll?ll L 2 (0,t;ff 3 (ng;R))' 

Thus, we are lead to 

\\(w,m% t <(c K +c s ti)\\(w,m% t +m^m% t 

+ Cs[N(u 0 ,(w i ) 3 i=1 ) + N((q i ) 2 =0 ) + M(f, Kg,nh)], 

which leads as in Section[§]to the introduction of a polynomial, this time of degree 
4, which does not bring any substantial change with respect to Sectional Note 
that the addition of C'«:||('U),q')||| t does not create any difficulty since a small K\ 
is chosen at the same stage as t \, and the conclusion is similar as in Section [9] 
from the continuity of ||(w;, g)||^ on [0, T„] which is established in the same way 
as the continuity of || {w, q)\\z t ■ We then infer that there is a time of existence of 
k for our smoothed problems, with a bound on || (u), g)||^ and thus on ||(ui, g)||z T 
independent of n. Existence follows then by weak convergence in Y? and uniqueness 
can be established similarly as for the compressible case in Section ITT! 
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